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I.  INTRODUCTION 


The  use  of  composite  materials  in  structural  design  has  been 
practiced  since  the  early  stages  of  the  human  civilization.  It  is 
only  lately  that  a  very  major  effort  has  been  made  in  improving  and 
comprehending  the  response  and  behavior  of  this  class  of  "tailorable" 
materials.  Depending  on  one's  needs,  modern  composites  may  consist 
of  a  number  of  different  types  of  constituents.  The  particular  class 
of  modern  composites  which  is  of  interest  here  is  the  one  consisting 
of  continuous  fibers  embedded  in  a  polymeric  matrix  or  binder. 
Structures  made  from  these  types  of  constituents  will  exhibit  some  form 
of  time-dependent  behavior  [1]  as  one  phase  (the  matrix)  is  time- 
dependent  and  the  other  (the  fiber)  may  or  may  not  be  time-dependent, 
depending  on  the  type  of  fiber  used  [2].  For  example,  automobile  tires 
reinforced  with  polyester  cord  will  correspond  to  the  former  and  epoxy 
based  composites  reinforced  with  graphite  fibers  is  representative  of 
the  latter.  In  addition,  irrespective  of  the  type  of  reinforcement, 
the  interface  between  the  fiber  and  matrix  is  time-dependent  £3].  For 
definitive  purposes,  only  epoxy  based  continuous  fiber  reinforced 
composites  are  investigated  here;  namely,  graphite/epoxy  laminates. 
Unless  otherwise  stated,  this  class  of  composites  is  the  one  addressed 
throughout  this  manuscript. 

When  considering  the  time-dependent  behavior  of  a  composite 
material,  there  are  at  least  four  accelerating  factors  (temperature, 
stress,  vibration  and  moisture)  [4]  that  will  affect  the  viscoelastic 
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behavior  of  a  material.  The  type  of  time  analogy  (superposition) 
selected  is  solely  determined  by  the  detrimenUl  effects  of  the  indi¬ 
vidual  accelerator  on  a  structure.  In  reality,  all  four  accelerators 
are  generally  present  in  most  structures.  Sirolation  of  experiments 
containing  all  four  accelerating  factors  and  boundary  value  definition 
of  these  experiments  are  indeed  cor^lex.  Consequently,  investigations 
pertaining  to  such  real-time  simulation  have  not  heretofore  been 
accomplished.  However,  individual  time  analogies  have  been  reported 
[5-10].  Unfortunately,  the  majority  of  these  limited  investigations 
are  neither  detailed  nor  complete  enough  to  be  comprehensively  applied 
to  two-phase  composite  systems  in  general.  The  primary  goal  of  this 
Investigation  is  to  address  these  latter  issues  as  outlined  in  greater 

detail  subsequently. 

The  objective  of  the  present  investigation  is  to  provide  an  in- 
depth  study  on  the  time- temperature  behavior  of  epoxy  based  continuous 
and  elastic  fiber  reinforced  composite  materials  and  explore  the  possi¬ 
bilities  of  using  the  time-terperature  analogy  as  an  accelerated 
characterization  method  to  predict  long-term  behavior.  Both  analytical 
and  experimental  time-temperature  procedures  will  be  presented.  The 
reason  for  considering  only  one  time  analogy  is  primarily  due  to  the 
constraints  imposed  by  time  limitations  which  would  rule  out  such  a 
diverse  investigation  as  multi -parametric  predictions.  However,  it  is 
felt  that  by  using  the  present  limited  study  other  environmental 
parameters  can  be  appropriately  and  confidently  incorporated  in  future 
multi -parametric  investigations  as  time  pemits.  The  reason  for 
selecting  the  time-temperature  analogy  is  that  among  the  four  analogies. 
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this  is  the  most  comnoo  and  also  the  most  detrimental  environmental 
parameter.  In  addition,  this  parameter  Is  one  to  which  epoxy  based 

composite  structures  are  always  subjected. 

Preliminary  investigations  on  the  time-temperature  behavior  of 
graphlte/epoxy  laminates  indicated  that  the  investigation  could 
essentially  be  separated  into  two  parts;  i.e..  room-temperature  and 
elevated  temperature  responses.  At  elevaud  temperatures,  especially 
at  or  above  the  glass-transition  temperature  (T,),  major  time-dependent 
behavior  was  observed.  At  room-temperature,  only  very  minor  visco¬ 
elastic  effects  were  observed.  Thus,  based  on  the  following  reasons, 
the  material  can  be  assumed  to  be  time  dependent  at  room  temperature; 

a)  In  the  analyses  of  Sandhu  [11],  Sendeckyj  et  al.  [12]  and 
others  [13],  where  only  quasi-elastic  behavior  is  assumed, 
good  correlations  between  analysis  and  experiments  were  ob- 

talned. 

b)  The  strain  rate  test  data  of  Daniel  et  al.  [14],  as  well  as 
VPI  data,  revealed  only  minor  viscoelastic  behavior.  How¬ 
ever,  the  results  were  inconclusive  due  to  the  scatter  in 

the  data . 

c)  The  creep  test  daU  presented  herein  shows  only  a  very  slight 
creep  behavior  for  off-axis  tensile  tests. 

For  these  reasons,  the  behavior  at  room  temperature  can  be  con¬ 
sidered  as  a  special  case  in  our  time-temperature  investigation.  For 
the  sake  of  completeness,  validation  of  test  methods  used  in  this  in¬ 
vestigation  and  to  provide  a  better  insight  to  the  present  and  future 
time-dependent  investigations,  this  special  case  is  summarized  and 
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substantiated  by  the  analytical  and  experimental  results  given  in 
Appendices  A  and  B,  where  a  synopsis  of  the  unnotched  behavior  is 
presented  in  Appendix  A  and  a  synopsis  of  the  notched  behavior  is 
presented  in  Appendix  B. 

The  fonnat  adopted  in  this  manuscript  is  to  first  present  the 
analytical  methodology  (Chapter  11).  In  this  chapter,  the  reduction 
of  the  anisotropic  time-temperature  constitutive  stress-strain  relations 
to  a  simpler  form,  a  discussion  on  the  time-temperature  superposition 
principle  and  the  creep  strength  prediction  methodology  are  presented. 
The  following  chapter  (Chapter  III)  presents  the  experimental  pro¬ 
cedure  including  our  observations  of  and  rationale  for  the  test  methods 
used.  In  Chapter  IV.  data  reduction  and  predictions  of  the  uni¬ 
directional  master  curves  obtained  from  short  term  tests  are  discussed, 
in  addition,  the  procedure  used  in  the  creep  to  rupture  predictions  is 
also  discussed.  The  analytical  and  experimental  results  are  compared 
and  discussed  in  Chapter  V.  Finally,  in  Chapter  VI  a  sumnary  of  the 
investigation,  as  well  as  conclusions  and  possible  avenues  of  future 
research  in  this  area  of  composite  materials,  is  discussed. 


II.  ANALYnCAL  METHODOLOGY 


The  vtscoelastic  structural  response  of  a  therniorheologically 
simple  material  (TSM)  is  affected  by  temperature  mainly  through  thermal 
expansion  and  changes  in  relaxation  and  retardation  times  (sometimes 
known  as  rheological  coefficients)  [15,16].  Thermodynamically,  the 
linear  viscoelastic  relations  bet^  stresses  and  strains  for  an 
anisotropic  TSM  can  be  written  compactly  by  means  of  the  Boltzmann 
superposition  integral  [17]  as 


_  a  (t  -  c') 


(1) 


and  in  the  inverse  form  as 


c,j  •  -  L  -  '•>  It 


(2) 


where 


o  and  7,  -  =  stress  and  strain,  respectively,  with  respect  to 
IJ 


’ij 


tine. 


coivliances,  respectively. 

?^j(C  -  f)  and  -  i’)  =  thermal  stress  and  strain 

characteristics  of  the  material,  respectively. 

aT  •  change  in  temperature.  T. 

r  ■  [  —  »  reduced  time 

)-.*T 


^ «  j  ^  »  reduced  time 
.  a.j(T)  »  shift  factor 
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and  i  ,j  »l^»l  "  ^  t2,3, - »6 

The  first  aed  second  lotesral  on  the  right  hand  side  of  eguatlons  (1) 
and  (2)  represents  the  response  due  to  mechanical  and  themal  loads, 
respectively,  and  the  lo»er  limit  (-)  Is  used  to  account  for  the  load- 
Ing  history  of  the  material.  Consequently,  equations  (1)  and  (2)  are 

sometimes  known  as  the  "Hereditary  Integrals." 

Under  isothermal  conditions,  equations  (1)  and  (2),  reduce  to 

-  ,  r  (r  -  c’)  — ^  dC 


respectively.  The  significance  of  the  material  properties  in  equa¬ 
tions  (3)  and  (4)  is  revealed  when  step-function  loads  are  applied. 
Mathematically,  the  step-function  stresses  can  be  represented  as 


where  are  constants  and 

0  for  c  <  0 

H(C)  -  •  „ 

1  for  c  >  0 

Substituting  aquations  (5)  and  (6)  Into  aquation  (4).  which  raprasants 
a  craep  tast.  and  for  virgin  matarials.  tha  resulting  equation  can  ha 


expressed  as 


^1j  *  ^ijkt^‘’^°kt 


The  inverse  of  equation  (7).  relaxation  test,  is 
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°fj  * 

(8) 

where 

!■  are  corsunts.  In  terns  of  the  real-time  parameter, 
k£ 

t,  equa- 

tions 

(7)  and  (8)  can  be  expressed  as 

(9) 

and 

(10) 

respectively.  At  this  stage,  either  equations  (7)  and  (8)  or  equations 
(9)  and  (10)  are  sufficient  to  represent  the  most  anisotropic  time- 
dependent  response  of  a  TSM  subjected  to  the  step-function  loads  des¬ 
cribed  earlier.  For  monolithic  materials  such  as  unidirectional 
composite  laminates,  the  thirty-six  material  constants  of  equations  (7) 
to  (10)  can  be  reduced  to  four  -hen  material  and  geometric  symnetries 
and  a  plane  state  of  stress  are  asst«d.  Thus,  equation  (9)  can  be 
expressed  as 

•‘22  *  ^26<‘>  “22’ 

Tij)  [5,6<‘>  ^26<‘> 

where  y^^2  “  ^^12*  simply  as 

and  the  inverse  of  equation  (12)  is 


ia)  •  CC(t)]{e*} 


(13) 
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It  should  be  mentioned  that  equations  (12)  and  (13)  are  expressed  in 
terms  of  the  global  coordinate  system.  If  the  local  and  global 
coordinate  systems  coincide  then  the  values  of  the  coupling  terms 
(S,g(t).  S2g(t),  C^glt)  and  ^gCt))  are  zero.  Tensorally.  the  local 
and  global  coordinate  systems  are  related  by  a  transformation  matrix, 
[TM],  which  is  expressed  as 


[TM] 


Cos^e  Sin^e 

7  2 

Sin  e  Cos  e 

Sine  Cose  -Sine  Cose 


-2Cos9  Sine 
2Cose  Sine 
Cos^e  -  Sin^e 


(14) 


where  e  is  the  angle  between  the  local  (1-2  coordinates)  and  the  global 
(x-y  coordinates)  axes,  as  shown  in  Figure  1.  The  relationship  between 
the  compliance  expressed  in  terms  of  the  local  coordinates,  [S(t)],  and 
the  compliance  expressed  in  terms  of  the  global  coordinates,  [?(t)],  is 

found  to  be 

[5(0]  «  [TM]*’  [S(t)][TM]  05) 


where 


[S(0]  = 


S<|^(t)  Si2(f) 
Si2(t)  S22(0 
0  0 


0 

0 


(16) 


where  S^^(t)  and  S22(0  are  the  compliances  along  and  transverse  to  the 
fibers,  respectively,  6^2^)  *  v.|2(t)S^l(t),  ^^20)  "is  the  major 
Poisson’s  ratio,  and  Sgg(t)  is  the  intralamina  shear  compliance.  Thus, 
knowing  these  four  tine-dependent  constants,  the  creep  properties  of 
any  arbitrary  unidirectional  laminate  can  be  obtained  from  equation  (15). 
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Among  the  nine  tenns  of  the  [5(0]  matrix,  the  first  diagonal  tenn 
{S^^(t).  axial  compliance)  is  pert\aps  the  most  important.  By  substi¬ 
tuting  equations  (14)  and  (16)  into  equation  (15)»  the  resulting  axial 
compliance  tem  is  expressed  as 

S„(t)  •  S,  (t)  ■  S„(t)  Cos^  *  (Sjj(t)  »  2S,j(t))  Cos^e  Sin^e 

.  (17) 


Equation  (17)  is  the  time-dependent  analogue  of  the  elastic  orthotropic 
transformation  equation.  In  terns  of  stresses  and  strains,  the 
principal  material  compliances  of  equation  (17)  are  expressed  as 


c,i(t) 

^1^^^  *  a\y 

(18a) 

II 

(18b) 

SfigCt)  »  44— 

66  x,2 

(18c) 

and 

-  li, 

(18d) 

Experimentally,  the  stresses  and  strains  of  equations  (18a)  and  (18d) 
are  obtained  from  a  [0“]^  uniaxial  test,  those  of  equation  (18b). are 
obtained  from  a  [90^^  uniaxial  test,  and  those  of  equation  (18c)  can 
be  obtained  from  a  [lO®]^  uniaxial  test  or  other  shear  testing  as 
described  in  Appendix  A. 

By  the  same  procedure,  the  other  five  elements  of  the  symmetric 
[?(t)]  matrix  can  be  obtained.  Thus,  by  the  application  of  appropriate 
assumptions  and  tensoral  arguments,  one  can  obtain  the  required 
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constitutive  eooetions  for  unidirectional  composites.  However, 
physically  the  time-temperature  verification  and  data  reduction  of  the 
constitutive  eouation  is  no  trivial  matter.  For  continuity  sake,  this 
aspect  will  be  placed  in  perspective  and  appropriately  discussed  in 
the  data  reduction  section  of  this  manuscript. 


Time-Tesperature-Superposition  Principle  (TTSP) 

Basically,  the  HSP  or  method  of  reduced  variables  [18,19]  is 
a  method  by  which  long-tem  time-dependence  can  be  determined  by 
short-tem  response  taken  at  various  elevated  temperatures.  It  is 

based  on  the  relation 


PawL 


(19) 


where  E,,(t.T)  and  E„(f,T,)  are  the  ..dull  at  temperature  T  and 
reference  temperature  T„,  respectively,  p  and  p„  are  the  densities  at 
temperature  ^  and  the  reference  temperature,  T^,  respectively.  For 
solids  the  ratio  p„/p  is  very  small  and  is  often  neglected.  The  time 
p.ra.«ters  f  and  t  are  related  by  the  shift  factor,  a^.  and  is  defined 

by  the  relation 

-ft 

It  must  be  mentioned  that  eouation  (19)  as  written  is  only  valid  for 
rela«at1on  moduli.  For  coimiosftes.  the  axial  creep  compliance  relation- 

ship  can  be  expressed  as 


Sxx<‘‘‘V‘Tr5)a<‘’''' 


(21) 


12 


or  in  the  reciprocal  form  as 


(22) 


Due  to  the  very  smll  change  In  densities  with  respect  to  tenperature. 
the  density  ratio  is  neglected  in  equations  (21)  and  (22). 

The  implication  of  equations  (19)  and  (22)  is  that  the  TTSP  will 
enable  one  to  reduce  the  complei  tine-teinierature  dependence  of  a 
viscoelastic  material  over  several  temperatures  to  a  single  time 
dependence  at  a  certain  reference  tanperature  and  a  separate  tempera- 
tore  dependence.  The  latter  dependence  is  represented  by  the  variation 
of  the  shift  factors  with  respect  to  the  isothermal  ambient  temperature. 
Physically,  it  would  nean  that  by  conducting  short-term  tests  at  dif¬ 
ferent  ambient  temperatures,  the  long-tenn  response  of  the  viscoelastic 
material  is  represented  by  a  corposite  curve  constructed  by  shifting 
all  of  the  short-tern  curves  to  obUin  a  master  curve.  This  master 
curve  is  obtained  by  the  horizontal  shifting  of  the  reduced  viscoelastic 
functions  with  respect  to  a  reference  temperature  along  the  logarithmic 

time  axis. 

As  such,  the  role  of  the  shift  factor  is  of  great  importance 
when  the  TTSP  is  invoked.  Generally,  for  artient  temperatures  below 
the  glass  transition  ten»rature  (T,).  the  shift  factors  are  related 
to  the  ambient  tenperature  by  an  apparent  activation  energy  (aH,) 
defined  as:  [18] 


d(Ln  ay) 

^a*"  aciTTF 


(23) 
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^ere  R  ts  the  eniversel  sat  constant,  U  Is  the  natural  logarithm 
and  the  other  parameters  are  as  previously  defined.  If  Ln  a^  is 
directly  proportional  to  1/T.  equation  (23)  can  be  further  simplified 
to  an  Arrhenius  type  equation,  that  is, 

ah,  •  2.303R  Logjp  a^  ]}  - 

The  implication  of  equation  (24)  is  that  the  activation  energy  Is 
constant  for  ambient  temperatures  below  the  T^.  If  the  ambient 
temperature  Is  at  or  above  the  T,  then  the  shift  factor  Is  related  to 
the  ambient  temperature  by  an  empirical  equation  known  as  the  Uilllam. 
tandel-Ferry  (ULF)  equation  [18-20],  The  WIF  expression  is  given  as 


where  and  Cg  are  material  constants. 

Senerally.  the  applicability  of  the  nSP  needs  no  verification. 
However,  due  to  morphological  changes  associated  with  changes  in 
temperature,  load  and  duration  of  the  test  [21,22],  poor  correlations 
between  the  deformation  predicted  by  nSP  and  experiment  may  occur. 
Often,  these  deviations  are  minimlied  by  shifting  the  deformation 
versus  log  time  curves  vertically.  In  addition  to  the  horizontal  shifts 
[23-26].  When  the  latter  procedure  Is  applied  to  the  construction  of 
the  master  curve,  the  material  Is  generally  considered  to  be  thermo- 
rheologlcally  complex  (TCH).  The  only  verification  of  the  applicability 
of  the  nSP  to  any  particular  material  Is  made  by  performing  long-term 


experiments. 
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Tine-Dependent  Strength  Prediction 
The  two  previous  sectloos  of  this  chapter  were  devoted  mainly  to 
the  short-  and  lono-term  time-temperature  behavior  of  unnotched 
composite  materials.  To  completely  characterize  the  time-temperature 
stress-strain  behavior  of  a  material  requires  that  a  criterion  be  used 
to  define  time-dependent  failure  or  more  specifically  time-dependent 
separation  or  rupture.  Analytically,  failure  criterion  or  criteria 
„st  be  invoiced  if  this  aspect  of  the  material  behavior  is  to  be  con- 

sidered. 

To  date,  a  larqe  nmnber  of  time-independent  failure  criteria  for 
anisotropic  materials  have  been  proposed  [27]  but  to  the  author's 
knowledge,  no  time-dependent  anisotropic  failure  criteria  have  been 
proposed.  An  alternative  would  be  to  apply  the  correspondence 
principal  [28]  to  the  time  independent  criteria  and  derive  its  time- 
dependent  analogue.  However,  such  a  procedure  would  require  the 
mathematical  functional  represenution  of  the  time-dependent  material 
properties.  A  variety  of  empirical  and  nwerical  techniques  for 
acquiring  the  viscoelastic  functions  have  been  proposed  [29-30]. 
Depending  on  the  material  behavior  and  the  degree  of  sophistication  of 
the  mathematical  representation,  the  failure  criterion  derived  in  this 

manner  may  or  may  not  be  appropriate  [31]. 

In  addition  to  the  cmnpleiities.discussed  above,  the  results  of 
Grinman  and  Gol'ctnan  [32]  showed  that  the  isotheimal  failure  behavior 
of  isotropic  high  density  polyethylene  varies  from  a  quasi-brittle  (at 
medium  and  high  stress-levels)  to  a  brittle  (at  low  stress-levels) 
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fracture.  The  implication  is  that  the  shape  of  the  failure  surface 
varies  with  time,  that  is.  from  a  criterion  with  dependent  failure 
modes  (short-tenn)  to  a  criterion  with  independent  failure  modes  (long¬ 
term).  Consequently.  Grinman  and  GoTdman  used  two  different  criteria 
to  correlate  their  creep  to  rupture  experimental  results  with  analysis. 
In  addition,  these  authors  assumed  that  the  principal  strengths  are 
functions  of  time  and  isochronal  principal  strengths  are  used  for  their 

predictions . 

In  view  of  the  limitations  in  using  the  correspondence  principal 
to  derive  time-dependent  failure  criterion  and  the  failure  mode  varia¬ 
tion  with  time,  discussed  above,  the  philosophy  adopted  here  in  pre¬ 
dicting  delayed  creep  failure  is  to  utilize  short-term  test  results  to 
predict  long-term  failure.  The  rationale  and  procedure  used  in  the 
predictions  will  be  discussed  in  subsequent  paragraphs. 

Consider  a  creep  to  rupture  test.  Assuming  that  the  time- 


temperature  variations  of  the  principal  creep  compliance  master  curves 
are  identical  to  the  variations  of  their  respective  principal  strengths, 
the  following  relation  between  strengths  and  compliances  is  postulated. 


that  is. 


RItl 


S  (t.T) 
nm 


where  X^(t.T)  «  creep  strength  at  T*C 
niii 

X  (taT.T')  *  creep  strength  at  T'®C 
nm  I 

m  ■  1.2.6  (sunination  not  implied) 


(26) 


and 


T®C  >  T'^C 
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It  should  be  mentioned  that  the  assumptions  used,  =  X^(t)  and 
the  identical  variation  of  the  principal  strengths  and  compliances,  are 
similar  to  those  used  in  reference  [32]  and  [33],  respectively.  The 
implication  of  equation  (26)  is  that  if  the  principal  compliances  obey 
the  TTSP  (that  is,  if  master  curves  of  the  principal  material 
compliances  can  be  obtained)  then  the  time-dependent  failure  envelopes 
will  possess  the  same  shape  as  the  master  curves.  Perhaps  more  im- 
portantly.  equation  (26)  assumes  identical  shift  factors  for  creep 
strengths  and  compliances ,  respectively.  The  physical  interpretation 
of  equation  (26)  is  that  if  one  principal  creep  strength  is  known,  the 
other  principal  creep  strengths  (within  the  time-temperature  range  of 
its  respective  corpliance  master  curve)  can  be  calculated. 

Once  the  three  principal  strength  master  curves  at  a  particular 
temperature  are  known,  the  application  of  the  three  time-dependent 
strengths  to  failure  prediction  of  different  plane  states  of  stress 
is  perfomed  by  using  isochronal  principal  strengths  in  conjunction 
with  a  time-independent  failure  criterion  derived  from  a  tensoral  point 
of  view.  Based  on  the  brief  but  concise  survey  and  evaluation  of  time- 
independent  failure  criteria  in  Appendix  A.  the  criterion  chosen  here 
is  the  one  proposed  by  Puppo  and  Evensen  [34].  Besides  the  advantages 
discussed  in  Appendix  A,  the  preference  of  this  criterion  over  the 
others  is  primarily  based  on  convenience.  That  is,  the  variation  of 
the  failure  surface  with  time  can  be  appropriately  accounted  for  by  the 
interaction  factor  and  its  exponent.  The  criterion  proposed  by  Puppo 
and  Evensen  for  a  plane  state  of  stress  is  expressed  as 
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where  y  * 


Tn 
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and 


n  is  a  material  parameter.  If  isochronal 


strengths  are  used  then 


and 
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Again,  the  strengths  in  equations  (27a)  and  (27b)  are  along  the  material 
axes.  For  a  prescribed  plane  state  of  stress  due  to  creep  loading,  the 
applied  stresses  are  constant  but  the  principal  strengths  are  time- 
dependent.  As  time  progresses,  the  principal  time-dependent  strengths 
decrease.  Failure  is  said  to  occur  when  the  sum  of  the  terms  on  the 
left  hand  side  of  equation  (27a)  or  (27b)  equals  to  unity. 

Thus,  from  the  Boltzman  superposition  integral  expressed  in  equa¬ 
tion  (2),  it  is  shown  how  one  could  utilize  the  TTSP  and  isochronal 
strengths  to  predict  the  time-temperature  behavior  of  a  unidirectional 
composite  laminate  subjected  to  creep  loading.  The  validity  of  the 
assumptions  made,  can  only  be  verified  by  experiments.  Only  one 
composite  system  is  chosen  to  verify  the  methodology  described  above. 

It  is  felt  that  if  the  experimental  results  do  indeed  correlate  with 
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the  analytical  representation  then  the  methodology  can  be  extended  to 
Other  composite  systems. 


III.  EXPERIMENTAL  CONSIDERATIONS  AND  OBSERVATIONS 

As  mentioned  earlier,  graphite/epoxy  composite  laminates  were 
used  in  the  experiMnUl  progrw-  This  particular  composite  system 
was  selected  because  of  the  stability  of  the  graphite  fibers  under  dif¬ 
ferent  theraal  conditions,  that  is.  they  do  not  change  morphologically 
as  is  seen  in  some  glass  fibers.  Except  for  the  AS.3501  graphite/ 
epoxy  laminates  used  in  son*  of  the  tests  in  Appendices  A  and  B.  all 
the  experimenUl  results  presented  here  were  obtained  from  T300/934 
graphite/epoxy  laminates,  the  fiber  volume  ratio  of  all  cured  panels 
were  approximately  605  to  655. 

Since  only  the  time-tercerature  analogy  is  considered  here,  the 
effects  of  the  other  three  accelerating  factors  (moisture,  vibration 
and  stress)  were  eliminated  or  minimized  in  the  experiments  as  follows; 

The  first  factor,  roistare.  was  eliminated  by  heating  the  speci¬ 
mens  in  an  oven  at  120«C  for  a  few  hours,  weighing  the  specimens 
individually  and  repeating  the  procedure  until  the  weight  of  the  speci¬ 
mens  stabilizes.  Generally,  tne  subilization  of  the  specimen's  weight 
was  observed  after  the  second  thermal  cycle.  These  -dry"  specimens 

were  stored  in  a  desiccator  until  tested. 

The  second  factor,  vibration,  was  avoided  by  performing  relaxa¬ 
tion  or  creep  tests.  Due  to  testing  requirements,  creep  testing  was 
preferred  because  it  was  easier  to  maintain  constant  stress  than 
constant  strain.  Also,  creep  represents  a  stress  state  that  results 
in  a  catastrophic  failure  mode  more  often  than  relaxation. 
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Finally,  the  third  factor,  stress,  was  minimized  by  operating 
at  stress-levels  such  that  linear  behavior  would  occur.  For 
example,  consider  the  plot  of  the  creep  stress-level  versus  15-minute 
creep  strain  of  Figure  2.  The  initiation  of  the  non-linear  behavior 
likely  occurred  at  approximately  the  point  where  the  stress-strain 
plot  begins  to  deviate  from  the  straight  line  passing  through  the 
lower  stress-level  points.  As  expected,  the  non-linearity  always 
Initiated  at  a  higher  stress-level  as  the  isothermal  ambient  tempera¬ 
ture  was  lowered.  By  operating  at  stress-levels  much  lower  than  the 
non-linear  Initiation  stress-level  of  the  respective  laminate,  the 
material  behavior  could  be  assumed  to  be  linear  viscoelastic. 

Another  precautionary  measure  suggested  by  Shapery  [26]  was  also 
adopted  In  the  investigation.  That  is,  under  mechanical  creep  loading, 
if  the  initial  jump  in  strain,  c^,  was  equivalent  to  that  when  the 
specimen  was  unloaded  the  material  behavior  was  assumed  to  be  linear, 
as  illustrated  in  Figure  3. 

To  insure  that  meaningful  data  were  acquired  (that  Is,  the  creep 
response  measured  was  primarily  due  to  the  specimen  and  not  due  to  the 
adhesive  used  to  bond  the  end-tabs  to  the  specimen  or  a  combination 
of  both  responses),  the  gripping  of  the  specimen  onto  the  testing 
machine  was  not  done  with  conventional  wedge  grips.  Instead,  a  modified 
version  of  the  gripping  system  suggested  by  Chamis  and  Sinclair  [35] 
was  used.  The  modification  made  here  was  the  Introduction  of  a  pin 
through  the  mid-section  of  each  clamped  section  of  the  specimen.  The 
advantages  gained  from  this  gripping  system  were:  The  In-plane  bending 
of  the  tensile  specimen,  discussed  by  Pagano  and  Halpin  [36],  and 
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slippage  In  the  gripped  sections  .ere  eliminated.  It  should  be  noted 
that  slippage  may  not  create  as  large  an  error  in  the  creep  tests  as 
in  relaxation  tests  but  could  still  lead  to  a  very  severe  consequence 

if  not  considered. 

In  addition  to  the  above  considerations,  specimens  need  to  be 
conditioned  thermally  and  mechanically  [37],  prior  to  data  acquisition. 
This  requirement  is  due  to  the  fact  that  epoxy  based  composite 
materials  are  normally  cured  at  a  temperature  below  their  vitrifica¬ 
tion  temperature.  Consequently,  the  properties  of  such  materials  may 
vary  if  the  material  were  cycled  mechanically  and/or  thermally  [38]. 
Such  behavior  may  be  observed  in  some  photoelastic  material  where  the 
stress-optic  constants  vary  with  themal  cycles.  Based  on  these  find¬ 
ings,  some  of  the  laminates  used  here  were  thermally  and  mechanically 

conditioned. 

After  five  thermal  cycles,  both  the  mechanical  and  thermal 
properties  were  essentially  the  same  for  the  fifth  cycle  as  the  first. 
Figure  4  shows  the  expansional  behavior  of  the  [90»]q^  laminate  for  the 
first  three  thennal  Cycles.  Three  distinct  values  of  the  coefficients 
of  thennal  expansion  in  the  90®  (transverse  to  the  fibers)  and  45 
directions  were  observed.  The  averaged  thermal  strains  were  measured 
by  two  rectangular  rosettes  mounted  on  each  side  of  the  specimen  and 
in-situ  temperature  cocpensated  by  two  identical  rosettes  mounted  on  a 
material  (known  as  ’Flat  Glass®  by  Coming  Glass  of  Blacksburg)  with  a 
thennal  coefficient  of  expansion  of  approximately  0.5  ue/®C.  The  tri- 
linear  thennal  behavior  of  the  material  occurred  in  the  three 
temperature  ranges  indicated  in  Figure  4.  Subsequently,  two  distinct 
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glass  transition  temperatures  were  ascertained  from  this  study;  the 
secondary  glass  transition  temperature  at  60“C  and  the  primary  glass 

transition  temperature  at  ISO^C. 

Upon  conpletion  of  the  thermal  cycles,  the  specimens  were 

mechanically  conditioned  by  subjecting  the  specimens  to  creep  loads  of 
approximately  50%  of  the  room-temperature  ultimate  load.  After  15 
minutes,  the  specimens  were  unloaded  and  allowed  to  recover  to  approxi¬ 
mately  zero  strain.  The  process  was  repeated  five  times  and  the  speci¬ 
mens  were  further  subjected  to  fifty  low  frequency  cyclic  loads.  At  the 
end  of  the  thermal  and  mechanical  cycles,  the  properties  of  the  speci¬ 
mens  were  essentially  the  same  as  they  were  in  the  virgin  state.  This 
tended  to  indicate  that  thermal  and  mechanical  cycling  were  unnecessary 
for  our  material  and  that  the  specimens  could  be  used  repeatedly  for  at 
least  the  number  of  thermal  and  mechanical  cycles  mentioned. 

After  establishing  the  essential  experimental  procedure,  16- 
minute  uniaxial  creep  tests  followed  by  recovery  to  approximately  zero 

strain  were  performed  on  the  ^8s’  ^8s’  ^8s’  ^8s* 

[60*]q^.  [75“]q5  and  tensile  specimens,  at  different  isothermal 

tenperature  levels.  The  reason  for  choosing  a  time  of  16  minutes  was 
due  to  the  experimenUl  procedure  adopted.  That  is,  data  were 
acquired  1  minute  after  the  initiation  of  the  load  and  at  3-minute 
Intervals  thereafter.  A  typical  creep- recovery  curve  is  as  shown  in 
Figure  3.  The  [IO'Iqs  laminate  was  used  to  obtain  the  time-dependent 
intralamina  shear  stress-strain  response  as  discussed  in  Appendix  A. 

Upon  completion  of  the  short-term  tests.  25-hour  continuous  creep 
tests  were  performed  on  the  same  specimens  at  180»C.  Complete  recovery 
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of  these  medium-tenn  tests  were  not  obtained  because  of  the  length  of 
.time  required  by  the  recovery  process.  The  purpose  for  this  senes 
of  tests  was  to  validate  the  TTSP  as  discussed  earlier. 

Finally,  new  specimens  were  used  to  ascertain  the  time  frame  of 
the  creep  to  rupture  (delayed  failure)  tests.  At  least  two  specimens 
of  identical  orientation  but  subjected  to  different  stress-levels  were 
used.  In  all  the  creep  to  rupture  tests,  the  ambient  temperature  was 
maintained  at  ISO^C.  It  should  be  mentioned  that  delayed  failures  were 
not  observed  at  all  the  prescribed  stress-levels  of  the  respective 
laminates.  The  specimens  that  did  not  fail  in  72  hours  were  replaced 
by  another  identical  specimen,  subjected  to  higher  stress-levels. 

In  all  the  creep  tests,  only  the  specimen  displayed  no 

creep  behavior.  At  elevated  temperatures  up  to  210“C.  the  mechanical 
properties  were  essentially  the  same  as  the  room-temperature  values, 
in  addition,  the  major  Poisson’s  ratio  {v^2*  obtained  from  this  laminate) 
is  found  to  be  time-independent.  This  observation  clearly  substantiates 
one  of  the  earlier  assumptions;  that  is.  the  fiber  was  fundamentally 


elastic. 


IV.  DATA  REDUCTION 

EXPERIMENTAL  MASTER  CURVE  DETERMINATION 

The  data  reduction  of  the  short-term  creep  compliances  was  per- 
fonned  by  measuring  the  isochronal  strains  from  1  minute  to  16 
minutes  at  3-minute  increments  and  dividing  the  six  measured  strains 
by  the  applied  stress  at  a  particular  temperature  level.  This  procedure 
was  repeated  for  all  temperature  levels.  The  compliances  were  in  turn 
multiplied  by  the  factor  where  T  is  the  absolute  temperature  at 
which  the  test  was  performed  and  T^  =  T^  »  453«K.  The  values  obtained 
in  this  manner  are  frequently  referred  to  as  reduced  values  and  form 
the  basis  of  the  HSP  as  discussed  in  Chapter  II.  The  primary  glass 
transition  temperature  was  taken  as  the  reference  temperature  mainly 
for  two  reasons:  First,  this  temperature  is  unique  for  each  polymeric 
material  and.  second,  such  a  selection  simplifies  the  task  of  verify¬ 
ing  the  Arrhenius  and  WLF  equations  discussed  in  Chapter  II.  By 
plotting  these  reduced  compliances  versus  log^g^^^"®^  joining  the 
shifted  Isothermal  reduced  values  by  a  smooth  curve,  segments  of  a 
long  continuous  curve  can  be  visualized.  Examples  of  such  plots  are 
shown  in  Figures  5  to  13.  except  that  the  reciprocal  of  the  reduced 
compliances  are  used  Instead  of  moduli.  The  purpose  of  this  procedure 

will  be  discussed  in  the  next  chapter. 

By  shifting  horizontally  the  reduced  curves  obtained  at  tempera¬ 
tures  above  the  Tg  to  the  right  and  the  ones  obtained  at  temperatures 
below  the  Tg  to  the  left,  a  smooth  continuous  master  curve  is  obtained. 
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The  shift  factors  (a^d))  are  equivalent  to  the  amount  of  horizontal 
shift  of  the  respective  isothemal  reduced  curves.  Portions  of  the 
respective  master  curves  are  shown  In  Figures  5  to  13.  It  should  be 
mentioned  that  this  is  a  "trial  and  error"  graphical  smoothing  pro¬ 
cedure  in  which  errors  may  arise.  However,  the  principle  yields 
imporunt  Information  over  a  broad  time  scale  as  indicated  by  the 
complete  master  curves  of  the  three  principal  material  properties  shown 
in  Figures  14  to  16.  In  these  figures,  the  time  scale  of  material 
response  has  been  expanded  to  over  twenty-six  decades  of  time. 

Master  Curve  Prediction 

From  an  experimental  and  computational  standpoint,  the  conven- 
tlonal  procedure  used  In  obtaining  master  corves  Is  a  very  tedious  and 
cumbersome  one.  This  procedure,  as  It  Is.  would  have  to  be  repeated 
for  every  arbitrarily  oriented  laminate.  Obviously,  from  the  designer’s 
point  of  view  and  for  economic  reasons,  this  Is  a  very  serious  draw¬ 
back  to  the  TTSP.  Thus,  a  new  technique,  based  on  the  principal 
material  master  corves.  1s  proposed  to  overcome  the  disadvantages  in  the 
conventional  master  curve  data  reduction  procedure  and  is  discussed  in 

detail  subsequently. 

Consider  the  orthotropic  transformation  equation  (Chapter  II) 
that  is, 

Sxx(t)  -  S,i(t)  Cos%  +  (Sgg(t)  +  2S,2^t))  Cos^e  Sin^e 
+  S22(t)  Sin^e 


(17) 


Master  Curvi 
at  lao-c. 
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For  the  type  of  t.«-phase  ™teri,ls  considered  herein  end  iron  oor 
enperi^entei  observation  on  these  »teriais  as  «11  as  those  of 
others  [37.39],  the  S„(t)  and  Sj^Ct)  terns  of  equation  (17)  are  foun 
to  be  tin-e-independent.  The  isochnmal  values  of  these  tuo  te™s. 
tabulated  in  Table  1.  do  indicate  ti«  and  temperature  independence 
at  least  within  the  t»^rature  ran,*  considered.  Thus,  equat.on  (17) 

can  be  further  simplified  as 

S„(t)  •  S„CosS  a  (S^(t)  a  3S,2)Sin^  Cos^^e  *  S2j(t)Sin*0  (28) 

.^re  tin*  dependent  te™s  are  indicate,  by  (t)  and  time  independent 
tenns  are  not  so  indicated.  The  difficulty  in  using  equat.on  (28)  m 
conlunction  with  the  four  asister  curves  is  that  reduced  master  corves 
of  the  S„  and  v.^  for  a  particular  t«,«rature  would  have  no  ^an.ng. 
Thus,  the  actual  constant  values  of  each  have  to  be  used.  On  the  other 
hand,  the  S,,(t)  and  S3,(t)  terns  are  «duced  values.  In  short,  equa¬ 
tion  (28)  will  not  be  closed  of  te™s  obtained  by  consistent  pro¬ 
cedures  if  the  actual  values  of  S,,  and  S,,  and  the  reduced  values  of 
S„(t)  and  S„(t)  are  substituted  directly  therein.  Any  ctmputational 
orocess  performed  using  equatiw  (28)  would  have  to  include  a  procedure 
to  mate  all  the  terms  consistent  with  each  other.  The  procedure  adopted 
was  to  use  the  actual  values  of  all  the  tenns  in  equation  (28)  as  will 

be  described  subsequently. 

Before  proceeding  into  the  computational  procedure,  it  would  be 
appropriate  to  discuss  a  few  relevant  and  essential  topics  on  the  time- 
ten^erature  behavior  of  two-phase  materials  in  general.  Consider  equa- 

tion  (22),  that  is. 


TABLE  1.  VALUES  OF  THE  RESPECTIVE  1-MINUTE  COMPLIANCES. 
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0.0468 

(0.0679) 


-0.0149 

(-0.0216) 


-0.0143 

(-0.0207) 


100 


180 


200 


210 


0.0411 

(0.0596) 

0.0404 

(0.0586) 

0.0402 

(0.0583) 

0.0400 

(0.0580) 


-0.0124 

(-0.0180) 

-0.0132 

(-0.0191) 

-0.0140 

(-0.0203) 

-0.0136 

(-0.0197) 


-0.0148 

(-0.0215) 

-0.0144 

(-0.0209) 

-0.0147 

(-0.0213) 

-0.0157 

(-0.0228) 


(22) 
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The  physical  interpretation  of  this  relationship  is  that  the  time,  t, 
required  to  reach  a  particular  response  (reciprocal  of  compliance  in 
this  case)  at  temperature,  T.  is  equivalent  to  the  time.  f.  required 
to  reach  the  same  response  at  temperature.  T^.  By  integrating  equation 
(20),  these  two  time  parameters  can  be  related  by  the  shift  factor  as 


t  -  a^f 
or 

Log,ot  •  Log^oa^  +  Log^pf  (W) 

The  main  difficulty  in  applying  equation  (17)  or  (28)  in  conjunction 
with  the  TTSP  to  two-phase  composite  systems  with  time-dependent  rein¬ 
forcements  is  in  the  utilization  of  the  shift  factors.  That  is.  the 
matrix  and  the  fibers  may  possess  different  shift  factors.  The 
representation  and/or  interaction  of  the  individual  shift  factor  of 
these  components  would  require  further  experimental  and  analytical 
study.  However,  in  two-phase  systems  containing  elastic  reinforcement, 
Moehlenpah  et  al.  showed  that  the  shift  factors  are  independent  of  the 
fiber  orientation.  Our  present  study  has  verified  that  observation. 
The  subsequent  computational  procedure  is  based  on  this  observation 
and  the  equivalent  response  relation  of  equation  (22). 

Consider  the  case  where  a  [e"]^  laminate  master  curve  reduced  at 
t'*C  is  required.  Since  the  temperature  dependence  is  represented  by 
the  variation  of  the  shift  factor  with  temperature  and  the  time 
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dependence  is  represented  by  the  master  curves  reduced  at  a  particular 
temperature.  T^.  the  physical  values  of  the  two  time-dependent  terms 
(S22(t)  and  SggCt)  of  equation  (28))  are  obtained  from  the  respective 
curves  and  the  required  master  curves  produced  in  the  following  manner: 

1.  For  the  particular  isothermal  ambient  temperature,  T^'C, 
a  generic  point  on  the  1/S22(t)  and  1/Sgg(t)  master  curves 
(reduced  at  T^'C)  must  be  established  first.  This  is 
accomplished  by  determining  the  shift  factor,  say  t,  for 
Tp*C  on  the  log^g^T  temperature,  as  indicated  in 

Figure  17a.  Using  this  shift  factor,  the  generic  point.  P, 
on  the  1/S22(t)  and  1/Sgg(t)  master  curves  are  located,  as 
indicated  in  Figures  17b  and  17c,  respectively. 

2.  The  two  reduced  time-dependent  prpperties  are  obtained  by 
measuring  the  ordinates  of  the  point  P  on  the  respective 

curves. 

3.  Obtain  the  physical  values  of  1/S22(t)  and  1/Sgg(t)  by 
using 

where  k  ■  2,6  (summation  is  not  implied).  It  should  be 
mentioned  that  T  »  t’  for  the  generic  point  and  is  equal  to 
T*  for  all  subsequent  time  increments,  as  indicated  in 

step  8. 

4.  Substitute  the  values  in  step  3  and  the  other  terms  of  the 
orthotropic  transformation  equation  into  equation  (28)  and 
compute  the  value  of 


Fiqure  17.  Schematic  for  Master  Curve  Prediction  Procedure. 


5.  Reduce  the  value  of  S^^(t),  obtained  in  step  4,  that  is, 

1  t; 

r-rb;;T“TT7^ 

XX  0  XX 

6.  the  next  valee  at  t  ♦  at,  tay.  of  the  time-dependent  proper¬ 
ties  are  obtained  by  using  t  ♦  at  to  locate  the  1/Sj2(t  *  at) 
and  l/S-£(t  ♦  at)  ordinates  In  Figures  17b  and  17c, 

DO 

respectively. 

7.  Use  the  t  +  At  value  to  obtain  the  equivalent  test  tempera¬ 
ture,  T’,  from  the  temperature  versus  log^g^T 

17d. 

8.  Repeat  steps  3  to  5  except  use  T  •  f  (In  steps  3  and  5)  and 
t  «  t  +  At. 

g.  Repeat  steps  6  to  8  for  every  additional  tine  Increment. 

By  plotting  the  reduced  1/S„(t)  values  obtained  In  the  above 
„„ner  versus  1og,„(t1me),  a  curve  representing  the  master  curve  of 
the  ta'lj  relo"'!  is  obulned.  this  procedure  can  be  used  to 

generate'any  arbitrarily  oriented  unidirectional  laminate  reduced  master 
curve  for  any  length  of  time  within  the  time  scale  of  the  input  data, 
as  well  as  any  arbitrary  ambient  temperatures  within  the  experimental 
range  of  the  principal  property  master  curves.  The  whole  procedure 
can  be  further  simplified  by  fitting  all  the  curves  in  Figure  17  with 
a  polynomial  (known  as  POLYRES  In  the  SSP  library  of  the  IBM  370/158 
computer  system)  and  used  as  input  data  for  the  computer  program 
written  to  perform  the  computations  described. 
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Strength  Predictions 

Uhen  predicting  the  delayed  failure  of  a  structure  subjected  to 
(say)  creep  loading,  there  are  essentially  two  interrelated  ways  that 
one  can  pose  the  probler..  One  way  is  to  determine  the  time  required 
for  the  structure  to  fail  when  it  is  subjected  to  a  known  stress-level. 

An  alternative  is  to  determine  the  stress-level  required  to  fail  the 
structure  at  a  given  tine.  Although  the  objectives  of  both  well  posed 
problems  are  similar,  the  computational  procedures  are  quite  different. 
Thus,  depending  on  one's  preference,  the  choice  is  solely  determined 
by  the  investigator.  Both  computational  procedures  and  an  alternative 
procedure  are  discussed  in  the  subsequent  paragraphs. 

Consider  the  case  of  determining  the  stress-level  required  to 
fail  a  specimen  at  a  particular  ambient  temperature.  T.  and  a  given 
length  of  tine.  t.  Sequentially,  using  the  analysis  described  in 

Chapter  II,  the  computations  are  performed  as  follows: 

1.  using  the  known  principal  strength  of  the  respective  laminates 

(obtained  from  tests  performed  at  T'),  calculate  the 
respective  isochronal  strengths  expressed  in  equation  (26). 

2.  Substitute  these  calculated  isochronal  strengths  at  T, 

into  equations  (27a)  and  (27b).  By  expressing  the  prescribed 
state  of  creep  stress  in  terms  of  the  local  coordinates  and 
increasing  the  applied  stress  incrementally,  the  failure 
stress-level  is  determined  when  the  sum  of  the  terms  on  the 
I.ft  hand  side  of  equations  (27a)  and  (27b)  equals  to  unity. 
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one  may  have  to  oerfom,  some  Iterations  (in  terns  of  the  applied 
stresses)  in  step  2,  in  order  to  achieve  unity.  For  uniaxially  applied 
stress.  the  incremental  procedure  of  step  2  can  be  avoided  and 
can  be  calculated  directly  as  follows:  Transform  o,  to  the  local 

coord 1 nates f 


•’ll  •  »,  » 

•’22  ■  ’x 

x-^2  *  ®x 


(30) 


Ho.  substitute  eouations  (30)  into  «iuations  (27a)  and  (27b)  to  obtain 
the  following  equations, 
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in  which  the  terns  are  as  previously  defined.  Thus,  using  these  latter 
equations  one  can  detemine  directly  the  required  uniaxial  stress-level. 

In  the  case  when  the  stresses  are  known  and  the  time-frame  is 
required,  increments  of  time  would  have  to  be  used  regardless  of  the 
types  of  in-plane  applied  stresses.  Except  for  sene  modifications, 
essentially  the  same  procedure  described  in  the  previous  paragraph  can 
be  used  for  this  prediction.  That  is.  the  isochronal  strengths  of 
step  1  is  required.  The  main  difference  is  in  step  2  where  the  applied 
stresses  are  constant.  By  taking  time  increments  and  substituting  the 
time-dependent  strengths  into  equations  (27a)  and  (27b).  the  required 
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time  frame  is  obtained  when  the  sum  of  the  terns  on 


the  left  hand  side 


of  the  failure  criterion  equals  to  unity, 
be  required  in  the  computations. 


Again,  some  iterations  may 


The  third  procedure  for  determining  either  the  time-frame  or 
stress-level  is  to  compute  creep  strength  master  curves  for  the  given 
sutes  of  stress  at  a  reference  temperature  using  the  procedure  des¬ 
cribed  for  tl«  first  case.  The  required  quantity  is  obtained  from  the 
time-dependent  failure  surface  directly  by  measuring  it  from  the  plot 


of  the  failure  surface.  The  main  disadvantage  in  the  simultaneous 
acquisition  of  both  quantities  is  that  different  surfaces  would  be 
«,„ir«i  for  different  ambient  temperatures.  However,  this  procedure 
is  i«ry  convenient  when  a  number  of  predictions  at  a  particular  ambient 


temperature  are  required. 


Y.  RESULTS  AND  DISCUSSION 


The  naster  curves  of  all  the  laminates  (except  the 
laminate)  obtained  from  the  short-tenn  (16-minute)  tests,  are  shown 
in  Figures  14  to  16  and  Figures  18  to  23.  The  results  in  these  figures 
represent  the  measured  master  curves  obtained  from  short-term,  16- 
minute  tests,  the  predicted  master  curves  obtained  using  the  computa¬ 
tional  procedure  described  in  the  previous  section  and  the  results  of 
our  25-hour  creep  tests  for  each  orientation  indicated.  In  all  these 
figures,  except  Figure  16,  a  linear  scale  was  used  to  represent  the 
magnitude  of  the  reciprocal  of  axial  compliance  instead  of  the 
logarithmic  scale  normally  used  in  the  master  curves  of  single-phase 
polymeric  materials.  This  is  primarily  due  to  the  small  volume  fraction 
of  the  epoxy  in  the  laminates  which  resulted  in  the  relatively  small 
variations  in  creep  compliances.  A  reason  for  using  the  reciprocal  of 
the  reduced  compliance  in  these  figures  rather  than  simply  the  reduced 
compliance  is  essentially  for  future  investigative  purposes.  That  is. 
a  number  of  authors  [1.37,39]  had  conjectured  that  for  temperatures 
below  the  Tg,  one  can  assume  the  behavior  of  the  material  to  be  quasi¬ 
elastic.  that  is 

Nomally.  in  single-phase  polymeric  materials  such  an  assumption  is 
invalid  [41].  The  only  means  of  validating  the  quasi-elastic  assump¬ 
tion  is  to  perform  relaxation  tests  on  the  same  material  and  compare 
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the  relaxation  results  with  the  reciprocal  of  reduced  compliance 
results  presented  here.  On  the  other  hand,  if  this  quasi-elastic 
assumption  is  valid  then  its  use  would  simplify  the  computation  of  the 
responses  of  general  laminates  considerably.  Primarily,  however,  we 
preferred  to  have  the  vertical  scale  be  comparable  to  moduli  rather 
than  compliance  simply  as  a  convenience  for  ourselves  and  for  others. 


As  indicated  earlier,  the  major  Poisson's  ratio  (v^2^  ''**5 
to  be  time- independent.  However,  the  minor  Poisson's  ratio  (v2p  ob¬ 
tained  fron  the  [SO’jgj  uniaxial  test)  was  time-dependent.  This  dif¬ 
ference  in  behavior  is  due  to  the  definition  of  Poisson's  ratio.  That 
is,  in  the  latter  case  it  is  defined  as 


“21 


CyyCt) 


The  experimental  results  of  the  laminate  indicated  that  the 

transverse  strain,  was  time- independent  and  the  axial  strain, 

e^(t),  was  time-dependent.  Thus.  V2^  was  found  to  be  a  time-dependent 
quantity.  Its  master  curve  is  shown  in  Figure  16.  Similar  results  for 
V2^  of  glass/epoxy  laminate  were  reported  by  Beckwith  [39].  However, 
he  attributed  such  a  result  to  experimental  scatter.  Based  on  the 
above  explanation  and  the  results  here,  perhaps  Beckwith's  conclusion 
was  incorrect. 

Except  for  Figure  16,  the  continuous  25-hour  tests  correlated 


well  with  the  predicted  master  curves  as  shown  in  Figures  18  to  23. 
Generally,  reasonable  agreement  was  obtained  for  all  the  three  results 
as  indicated  in  Figures  14  to  23.  The  large  deviations  in  the 
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Short-term  s'")  "suits  (Figures  21  and  23.  respective¬ 

ly)  may  have  been  due  to  experimeotal  and/or  the  graphical  procedure 
used  to  construct  the  master  curve  from  the  short-term  results.  A 
comparison  of  the  1/S„(1-minute)  results  obtained  from  the  transforma¬ 
tion  eguation  and  experiment.  Figure  24,  Indicates  that  for  these  two 
particular  laminates  deviations  occur  at  elevated  temperatures  of  about 
180'C  and  above.  The  large  deviations  of  the  tdS'lgj  and  C75”]gj 
master  curves  tends  to  substantiate  the  earlier  connent  on  the  con¬ 
ventional  procedure  used  in  constructing  master  curves  from  short-term 
tests,  that  is.  the  graphical  smoothing  procedure  may  create  errors. 
These  deviations  also  imply  that  the  master  curve  prediction  computa¬ 
tional  procedure  proposed  here  is  not  only  less  error  prone  but  is  a 
very  cost-effective  technique  for  the  accelerated  characterization  of 

composite  materials. 

Another  advantage  of  the  master  curve  prediction  procedure  is 
that  master  curves  of  off-axis  specimens  reduced  at  any  arbitrary 
reference  temperature  can  be  produced  at  a  relatively  small  cost  of  less 
than  a  minute  of  computer  time.  For  example,  consider  the  master 
curves  generated  tor  the  [SO'lgj  laminate  shown  in  Figure  25.  These 
master  curves  were  generated  for  the  three  reference  temperatures  of 
150*Cp  180®C  and  200®C  with  the  aid  of  the  shift  factor  versus 
temperature  plot. 

The  temperature  dependence  of  the  material  studied  here  is  shown 
by  the  Log,ga^  versus  temperature  plot  of  Figure  26.  The  deviations 
of  the  V2,(t)  results,  at  temperature-levels  less  than  the  T^,  from 
the  other  two  principal  properties  results  is  primarily  due  to  the 
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Figure  26.  Loq^q  a.|.  Versus  Temperature. 
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relatwel/  soli  transverse  strains  produced  by  the  uniaxially  stressed 

speci*”- 

tM-phase  raterials.  It  should  be  mentioned  that  the  ujift'  tesu'b* 
are  not  used  in  any  of  the  predictions.  The  reasons  for  presentin, 
results  for  vj,(t)  are  mainly  to  validate  the  elastic  symnetry  of  the 
aaterlal  (Table  1,  »here  $,2(1  minute)  =  S2)(l  minute))  and  for 
completeness  sake. 

Also  sho-n  in  Figure  26  is  the  correlation  between  the  HLF  equa¬ 
tion  and  the  experimental  results.  Good  correlations  were  obtained 
for  teaperatures  ranging  from  slightly  below  the  Tg  to  210»C.  The 
values  of  the  material  constants  used  to  obtain  this  correlation  were 
C  »  8.39  and  C2  *  52.97.  as  expressed  in  equation  (25).  Figure  27 
iLicates  that  log^Qa,  as  a  function  of  1/T  is  non-linear  for  tempera¬ 
tures  below  the  T^.  This  tends  to  invalidate  equation  (24)  which 
assimies  a  linear  function  and.  for  this  reason,  equation  (23)  would 
have  to  be  used  to  detemine  apparent  activation  energies. 

The  effect  of  terperature  on  the  strength  of  the  material  in¬ 
vestigated  here  is  shown  in  Figure  28.  The  experimental  strengths  of 
the  respective  laminates  were  obtained  by  ramp  loading  the  specimens  to 
failure  at  0.050  inch/minute  (0.127  cm/minute)  head-rate  and  at  the 
tMo  isotherwl  anbient  temperatures  indicated  in  the  figure.  This 
test  is  sometimes  assuned  to  be  a  very  short-term  creep  to  rupture  test 
[32].  The  predictions  were  made  using  the  failure  theory  of  Puppo- 
Evensen  which  was  discussed  earlier.  Reasonable  correlations  were  ob¬ 
tained  between  experiments  and  predictions  as  shown.  It  should  be 
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mentioned  that  by  osing  the  value  of  0.1  for  the  exponent  of  the 
Interaction  factor  of  equations  (27a)  and  (27b)  in  the  180'C  tests, 
a  slightly  better  correlation  was  obtained  as  opposed  to  using  the 
value  of  unity.  However,  the  value  of  0.1.  for  the  exponent  of  the 
Interaction  factor,  gave  a  much  better  correlation  with  the  delayed 

failure  test  results. 

The  third  data  reduction  procedure  for  predicting  delayed 
failures  and  a  value  of  0.1  for  the  exponent  of  the  interaction  factor 
Of  equation  (31)  were  used  to  produce  the  analytical  results  of  the 
respective 'laminates  shown  in  Figure  29.  together  with  experimental 
results.  The  open  symbols  represent  the  tests  where  delayed  failures 
were  observed  and  the  solid  symbols  represent  tests  where  the  speci¬ 
mens  did  not  fail.  The  solid  lines  represent  the  analytical  results. 

All  the  specimens  were  subjected  to  uniaxial  tensile  load.  The  three 
orthogonal  axes  represent  a„.  032.  and  time,  respectively.  The  ramp- 
loaded  to  failure  results  at  180“C  are  represented  by  the 
plane  at  10°  minute  and  are  assumed  to  be  the  very  short-term  creep  to 
rupture  tests.  All  the  open  symbols  either  coincide  with  their  respective 
strength  master  curves  or  are  in  the  vicinity  of  it.  The  solid  symbols 
are  always  below  their  respective  strength  master  curves.  In  addition, 
the  variations  of  the  failure  surface  with  time  are  compared  in  this 
figure.  That  is.  the  failure  surface  appears  to  "shrink"  as  the 
duration  of  the  tests  increases. 

A  better  representation  of  the  experimental  and  predicted  delayed 
failure  results  are  tabulated  in  Table  2.  In  this  table,  for  laminate 
orientations  equal  to  and  less  than  45“.  reasonable  correlations  were 
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TABLE  2.  COMPARISON  OF  EXPERIHEfiTAL  AND  PREDICTED  AXIAL 
CREEP  STRENGTH  AT  180*C. 


Duration. 
Laminate  of  Test 
Orientation  (ninutes) 


ns-iss 

C30«]o, 


[30-L 


C45*], 


tW]o, 


[60*], 


C75-], 


[75*3, 


586.0 


136.0 


613.0 


10A.O 


1230.0 


115.0 


2304.0 


Experimental 
X  ksi  (MPa) 

14.286 

(98.502) 

9.146 

(63.062) 

7.317 

(50.451) 

4.699 

(32.400) 

4.293 

(29.600) 

3.976 

(27.415) 

3.171 

(21.864) 

3.090 

(21.306) 

2.584 

(17.817) 


Predictions 
touation  (ilaj  Equation  13IPJ 


muabiwii  — ' - -  '  . 

X  ksi  (MPa)  *  ksi  (MPa) 


20.072 

(138.396) 

18.648 

(128.578) 

22.664 

(156.268) 

7.088 

(48.872) 

6.000 

(41.370) 

4.081 

(28.139) 

3.400 

(23.443) 

3.088 

(21.292) 

2.530 


14.971 

(103.225) 

7.811 

(53.857) 

5.513 

(38.012) 

4.058 

(27.980) 

3.623 

(24.981) 

3.114 

(21.471) 

2.589 

(17.851) 

2.632 

(18.148) 

2.108 
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obtatned  -ben  U>e  second  failure  criterion,  equation  (31b),  was  used. 

The  first  criterion,  equation  (31a).  qaye  excellent  correlations  for 
laminate  orientations  greater  than  45".  This  would  imply  that  for  the 
material  used  here,  the  first  criterion  would  be  appropriate  for 
e  >  45"  and  the  second  criterion  is  suitable  for  45”  i  o  s  0". 

Besides  the  good  overall  correlations  between  the  predictions 
and  experiments  on  the  delayed  failures.  Figure  29  clearly  illustrates 
the  dependence  of  the  creep  strength  master  curves  on  the  state  of 
stress.  That  is.  for  orienutions  between  0"  and  90".  the  location  of 
the  creep  strength  master  curve  on  the  envelope  is  deterained  by  the 
principal  stresses  (the  stresses  along  the  local  coordinates).  These 
principal  stresses  are  dependent  on  the  laminate  orienutions.  Conse¬ 
quently.  one  would  have  tt  produce  creep  strength  master  curves  for 
all  the  required  orienutions  either  from  experiments  or  other  means. 
From  the  designer's  point  of  view,  a  rational  and  convenient  methodology 
Should  be  available  u  predict  delayed  failures.  Based  on  the  results 
here,  it  is  felt  that  the  methodology  proposed  can  be  used  confidently 
and  iasily  to  generau  either  creep  strength  master  curves  of  any 
arbitrary  unidirectional  laminate  or  time-dependent  failure  envelopes 

at  the  required  ambient  temperature. 

Finally,  the  ertrittlement  of  the  maurial  due  to  the  duration 

of  the  ust  is  also  clearly  illustrated  in  Figure  29.  Consider  the 
straight  dashed  lines  connecting  points  A.  B  and  A'.  B'  as  indicated 
in  Figure  29.  For  very  short-term  creep  to  rupture  tests  (1  minute), 
the  material  failed  in  a  ductile  manner  as  indicated  by  the  degree  of 
curvature  of  the  surface  with  respect  to  line  AB.  Over  a  long  term 
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(10*  .inuttO.  O*  <le,ree  of  corvotore  of  tho  oorface  wUh  respact  to 
A-B'  decreases.  This  would  ia^.1y  a  n«ro  brittle  failure  at  10 
„i„utet  tbao  10“  .ioute.  This  type  behavior  Is  similar  to  the  ores  ob¬ 
served  for  sihsle  phase  material  subjected  to  different  biaxial 
sUtes  of  stress  132].  The  implication  is  that  methodologies  currently 
being  used  for  unfilled  polymeric  materials  can  be  extended  to  epoxy 

bas6d  composite  naterials. 
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VI.  CONCLUSIONS  AND  FUTURE  CONSIDERATIONS 

The  present  Investigation  has  been  concerned  with  the  time- 
tenperature  behavior  of  unidirectional  graphite/epoxy  laminates  at 

and  elevated  temperatures  while  general  laminates  have  been  only 
concerned  with  room  temperature  behavior.  For  elevated  temperatures. 

the  investigation  can  be  sumnarized  as  follows: 

.  The  expansional  behavior  of  the  material,  at  temperatures 
ranging  from  22»C  to  210%.  is  tri-linear.  Consequently, 
two  "knees"  are  observed  to  occur  at  60%  and  180%.  The 
first  "knee"  (at  60®C)  is  assumed  to  be  the  secondary  and 
the  other  "knee"  (at  180%)  is  assumed  to  be  the  primary  T^. 

.  The  viscoelastic  linearity  check  indicated  the  material 

behavior  to  be  linear  at  low  stress  levels  and  non-linear  at 
higher  stress-levels.  Our  material  was.  therefore,  thermo- 
rheological  ly  complex  but  a  thermorheologically  simple 
behavior  was  assumed  due  to  the  low  stress-levels  applied. 

.  Thental  and  mechanical  conditioning  of  the  material  used 

here  was  found  to  be  unnecessary. 

.  Among  the  four  principal  material  properties,  two  (S„  and 

time-independent  and  the  others  (^£2^^^ 

and  Sgg(t))  were  found  to  be  time-dependent. 

.  Master  curves  of  the  off-axis  specimens  were  obtained  by 

using  the  TTSP. 
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.  using  tl»  four  principal  material  properties  and  TTSP.  master 
curves  of  off-axis  specimens  were  produced  from  the  compute- 
tional  scheme  proposed  herein. 

.  The  predicted  master  curves  generally  agreed  better  with  the 
25-hour  continuous  test  than  the  ones  obtained  from  the  short- 

term  tests. 

.  Delayed  failures  -ere  observed  for  all  the  off-axis  specimens. 

.  A  cowoUtional  scheme  for  producing  master  creep  strength 
curves  of  any  arbitrarily  orientated  unidirectional  laminates 

at  any  reference  temperature  was  proposed. 

.  Heasonable  to  good  correlations  were  obtained  from  the  experi- 
nental  and  predicted  creep  to  rupture  results. 

.  The  master  creep  strength  curves  were  dependent  on  the  state 

of  stress. 

.  Depending  on  the  duration  of  the  creep  to  rupture  tests,  the 
„teHal  failure  mode  varied  from  ductile  (1  minute)  to  quasi- 
brtttle  (10*  minutes).  Both  experimental  and  analytical 
results  substantiated  the  embrittlement  behavior. 

Unnotched  and  notched  quasi-elastic  investigations  had  been 
teported  and  presented  in  the  appendices  including  a  sunmary  and  con- 

elusion  and  will  not  be  repeated  here. 

On  the  whole,  the  agreement  between  analyses  and  experiments  -ere 

satisfactory.  The  methodology  outlined  here  was  essentially  Intended 

to  be  used  for  the  accelerated  characterization  of  epoxy  based  con- 

timous  and  elastic  fiber  reinforced  composite  materials.  This  new 

«thodo1ogy  is  only  the  preliminary  portion  of  an  investigation  of  a 


major  effort  undertaken  for  the  purpose  of  evaluating  the  possibility 
of  forecasting  the  effects  of  the  four  accelerators  (vibration,  stress, 
moisture  and  temperature)  on  the  time-dependent  behavior  of  epoxy  based 
composite  materials.  However,  before  the  larger  issue  can  be 
addressed,  additional  verifications  and/or  modifications  on  the 
present  time-temperature  methodology  are  required.  Namely. 

1)  The  master  curves  presented  herein  were  all  reduced  at  180"C. 
Longer  tenn  continuous  testing  at  180%  and  a  few  other  lower 
isothermal  ambient  temperatures  should  be  performed.  The 
results  of  these  additional  tests  will  not  only  validate  the 
predictions  of  the  proposed  master  curve  computational  pro¬ 
cedure  but  also  the  applicability  of  HSP. 

2)  More  delayed  failure  tests  should  be  performed  at  lower 
ambient  temperatures  and  for  longer  time  periods.  The  results 
of  such  tests  would  be  useful  to  further  validate  the  delayed 
failure  prediction  methodology. 

3)  Relaxation  tests  should  be  done,  in  order  to  deteraine  whether 
the  quasi-elastic  assumption 

is  valid  at  temperatures  below  the  T^.  In  addition,  these 
tests  can  be  used  to  determine  the  interrelationship  between 
the  creep  compliance  and  relaxation  moduli  [42]  if  different 
than  given  in  the  above  equation.  Such  a  relationship  would 
in  turn  be  needed  to  determine  the  creep-relaxation  behavior 
of  cured  general  laminates. 
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4)  Develop  an  incremental  time-temperature  computational  method 
(similar  to  the  one  in  Appendix  A)  to  determine  the  behavior 
of  general  laminates.  The  methodology  developed  should  con¬ 
sider  creep  and  relaxation  responses  and  should  be  verified 
experirentally  along  the  lines  of  the  first  two  reconmenda- 

tions. 

5)  Explore  the  possibilities  of  a  more  effective  time-temperature 
failure  criterion. 

The  above  items,  when  completed,  would  subsequently  lead  to  a 
stress-dependent  HSP.  Before  attempting  to  use  such  a  joint 
•thodology  (referred  to  as  the  time-stress- temperature  superposition 
(TSTS)  [11]) »  the  time-stress  superposition  principle  (TSSP)  should  be 
verified.  At  present  there  are  at  least  two  approaches  [8,26]  on  the 
TSSP.  These  approaches  should  be  reviewed  and  experimentally  verified 
to  determine  their  feasibility.  The  methodology  developed  from  this 
study  should  be  incorporated  into  the  TTSP  and  the  resulting  methodology, 
TSTS,  should  be  verified  experimentally  using  the  rationale  and  proce¬ 
dure  described  in  the  manuscript.  Subsequently,  the  constitutive 
equations  used  in  the  TSTS  methodology  can  be  used  in  the  analysis  of 
the  notched  behavior  of  composites  subjected  to  different  stress  and/or 

temperature  levels. 

Finally,  the  confidence  and  experience  gained  from  the  TSTS 
.thodology  woTd  enable  one  to  incorporate  the  other  two  time  analogies 
(ti«-»ibration  and  time-moisture)  into  the  multi -parametric  predictions, 
again,  a  syrtiosis  between  analysis  and  experiment  must  exist  if  the 
resulting  methodology  is  to  be  meaningful. 


. 

VII.  BIBLIOGRAPHY 

1. 

H;tlDin  J.  C..  "Introduction  to  Viscoelasticity",  Corposite 

Materials  Hortshop,  Progress  in  Materials  Science  Senes,  Vol.  , 

1968. 

2. 

vauHoa  y  Y  '  Osipova,  E.  S.;  Faberova,  I.  !•*  snd  Artanova, 

K  Ya  ^"Tine-TefTperature  Superposition  Applied  to  the  Relaxation 
Lpertles  “  a  GUss-Reinfor«d  Plastic  and  Its  Katrix". 

Polymer  Hechanics,  March^April,  1972. 

I  ■  ■ 

3. 

Rov  A  .  and  Murthy,  P.  N.,  "Viscoelastic  Analysis  of  the  Residual 
Stms'in  ^ass  Fib^r-Reinforced  Plastics",  Fibre  Sncnce  and 
Technology,  Vol.  9,  1976. 

!  ■ 

4. 

Urthwctse,.  Tu.  S..  "Prediction 

of  Polyneric  Materials",  Polymer  Mechanics,  May-June,  19/2. 

5. 

u.inin  j  C  and  Paoano,  N.  J.,  "Observations  on  Linear 

Anisotropic  Viscoelasticity".  J.  of  Composite  Materials,  Vol.  2. 

1968. 

h 

1 

6. 

Urzhuratsev.  Yu.  S.,  and  Maksimov.  R.  0.,  "Multi -para^tric  Predic¬ 
tion  of  the  Creep  of  Polymer  Materials".  Polymer  Mechanics,  May- 

June,  1970. 

U*  ■ 

1'  •  , 

i. 

7. 

Turtsin'sh.  R.  P.,  Urzhumtsev,  Yu.  S..  and  Maksimov,  R.  0., 

•Static  and  Vibrational  Creep  of  a  Fabric-Reinforced  Plastic  in 

Shear  in  the  Plane  of  Reinforcement".  Polymer  Mechanics,  July- 
August.  1971. 

8. 

Urzhumtsev.  Yu.  S.,  "Time-Temperature  Superposition  for  Theiro- 
rheologically  Complex  Materials",  Polymer  Mechanics,  March-April, 

1974. 

9. 

Sims.  0.  F..  and  Halptn.  J.  C.. 

Elastic  and  Viscoelastic  Response  of 

Composite  Materials:  Testing  and  Design  (Third  Conference), 

ASTM  STP  546.  1974. 

10. 

Maksimov  R  0..  Mochalov.  V.  P..  and  Urzhumtsev,  Yu.  S.,  "Time- 
Moisture’superposition".  Polymer  Mechanics,  September-October, 

1972. 

11. 

Sandhu,  R.  S..  "Nonlinear  Behavior  of  Unidirectional  and  Angle- 
Ply  Laminates",  Journal  of  Aircraft.  February,  1976. 

74 

'  .  '  .  V  V.-  \  ‘  '  .  ■  '  * 

75 


«  r  **pf*3CtliT*C 

12.  Sendeckyj,  6.  P.»  ^^^Q5/5|o8’GrapMte/Epoxy*^Laminates--Part  1 : 

^ .  ..CM  "A  Method  of  Predicting  Non- 
liJiJr  B;ha;ior"of“Laminated  Composites",  J.  of  Composite 
Materials,  January,  1^69. 

«c4.Min  Rate  Effects  on  Mechanical 
SS!,^itiesSf’Frber''ion5osites",  HT  Research  Institute,  Report 
Mo.  IITRI  06073-lV,  June,  1976. 

Mechanics. 

July-August,  1968- 

j  .  CM  "Stress  Analysis  for  Linear  Visco- 
«•  rsti?-,itel!ia,r2i^trTe^e"raiu4riatioh’;.  Prans.  of  Soc.aty 
of  Rheology,  I960. 

”•  r“cK-a"i 

J.  Applied  Physics,  Vol.  35,  I9b4. 

a  1  D  -Viscoelastic  Propertios  Of  Poljnners-,  John  Hiley 

'®-  i'S.  ie«;d  Mfcioo,  New  York.  1970. 

.  V..  s  "Time-Temperature  Superposition  Review", 

Ii;5^r4jhlnic!;-January-February.  1975. 

>  S  a  V  -Properties  and  Structures  of  Polymers".  John 
*“•  KTso«.''i«.!™4:  Jork.  I960. 

a.  ,.rtene,.G.M..  and  Zelenev.  YU.  V..  Polymer  Mechanics.  HO.  . 

1971 

«.  i^lotko,.  A.  p..  zelenev.  Yu.  V..  and  Bartene..  6.  H..  Polisher 
Mechanics,  No.  6,  1971. 

,  .  «  r  "Time-Temperature  Superposition 

23.  Bahadur,  S.,  ?;;‘*5^“J®"?;-op;rtUs  of  the  Ethylene-Propylene 

^'rljs’^m^rr  Killed  Polymer  Science.  Vol.  16.  1972. 

’  March-June, 

1974.  .  _ 

^t^rlss^AnSogia  roisfrSrtfnrSni?Ierru;Ils"!'^ 

Mechanics,  May-June,  1974. 


76 


OC  «:rh;ir«.rv  R  A  ”0n  the  Characterization  of  Nonlinear  Visco- 
e)«t?7"atcriais".  Polymer  EnglneeP'PS 

27.  Seodhu.  «.  S.,  ■A.5Pr^»!/*J;“^'JS5^:^R!7°2^KT^‘’^o^’ce'' 
ni^r^murub!.  irlghtlpattersoo  A1r  Force  Base.  Ohio. 

28.  Flugge.  W..  •Viscoelasticity",  Sprioger-Verlag,  New  York.  1975. 

”•  ?„^12i.UstiJsfresrA:“ysfs"'1ro«edlogrIn.S?'c^^^^^^ 

for  Applied  Mechanics,  ASME,  1962. 

lir„ilrP^otovi;coe?as?rc"J;teHai ^•7'ExpeHienurMerhaiicl! 

SepteB<>er,  1964. 

31.  Wang,  A.  S.  0.,  Orexel  University.  [Private  Conwnicationl. 

“■  M^;iJM"inJl«iaVTMSto5\’pil’ym^™MeIh“ics.  toy-June, 

1976. 

F«cw™''p?SnH«  in’P0l^e“r"  ?IlidI:’'pSeSa°^d 
5:^7i;tt"c?eJ«  Jublishers.  John  Wiley  S  Sons.  Inc.  1964. 

£S“rian  Ler2o^rerst"res?;^.:«^^^^^^^^ 

S!Sr'p.Jierti.rof'nblVc'Jpis?tts"!°il!s“"  sJrin?K^^ 

Dallas.  Texas,  May  1977. 

“•  JJ^*T«tinj  rASlSl^pic  Liel"? rof  SLwtiw  Sner^U, 
January  1968. 

..  j  c..K.rva,.u  D  A  " VI scobI as ti c  Characterization  of 

”•  i°lliilin^;  ?;t,r-Relnfo«ed  Plistic".  0.  of  Composite  Materials, 
April  1971. 

,  u  DA  "^.irvAv  on  Time-Dependent  and  Fracture  Behavior", 

“•  SattH-alsrTes’unr^d  oSsign.  5th  Conference,  hew 

Orleans,  La.,  March  20-22,  1978. 

«f  .lANAAF.  San  Francisco,  February  26,  1976. 


K 


77 


An  ii««hio#tnAh  A  E  Ishai,  0.,  and  Di  Benedetto,  A.  T.,  ’The 
Effect  of  Temperature  on  the  Mechanical  Behavior  of  Epoxy 
Coeposites",  Polymer  Science  and  Engineering,  March,  19  1. 

41  ThMcaris  P  S.,  "Creep  and  Relaxation  Contraction  Ratio  of 
vucoiuaic  tta«naU".  J.  of  Mech.  Phys.  Solids.  »ol. 

12,  1964. 

47  De  Runtz  J.  A.,  Jr.,  and  Crossman,  F.  W.,  "Time  and  Temperature 
Effects  in  Laminated  Composites",  Proceedings  of  Computer 
Siojlation  of  Materials  Application,  Nuclear  Materials,  Vol.  20. 

1976. 

43.  Ashton.  J.  E..  and  Whitney.  J.  M..  "Theory  of  Laminated  Plates’. 
Progress  Kat.  Sc.  Series,  IV,  Technomic  Pub.  Co.,  1970. 

44.  Jones.  R.  "Mechanics  of  Composite  Materials",  McGraw-Hill 
Book  Company. 

45  Yeow.  Y.  T..  and  Brinson,  H.  ♦ ."Stress-Strain  and  Failure  ^ 

Properties  of  Graphite/Epoxy  Laminates  ,  VPI-E  76  , 

S.U..  Sept..  1976. 

46  Yeow  Y.  T..  and  Brinson.  H.  F.,  "An  Investigation  on  the  Tensile 
Moduli  and  Strengths  of  Graphite/Epoxy  Laminates  ,  Experimental 
Mechanics.  Nov.  1977. 

47  Yeow  Y  T..  and  Brinson,  H.  F..  "A  Comparison  of  Simple  Shear 
*  Characterization  Methods  for  Composite  Laminates  .  Composite, 

Jan.  1978. 

48.  RwlaHls.  R.  E..  -Flow  add  Fallura  of 

Inelastic  Behavior  of  Composites.  A.S.M.E.,  AMD-Vol  u.  iv/s. 

49  a  Ashkenazi.  E.  K.,  "Anisotropy  In  the  Strength  Co'^^truc^on 

Katerial",  C.  M.  Kirov  Wood  Tech.  Academy,  Leningrad,  May 

1961. 

b  Ashkenazi.  E.  K..  "The  Construction  of  Limiting  Surfaces  for 
Biax^  Stressed  Condition  of  Anisotropic  Materials  . 
Zavodskaya  Laboratorlyz.  Leningrad  Forestry  Academy.  Feb.  1964. 

c.  Ashkenazi.  E.  K.,  "Problems  of  the  Anistropy  Strength’. 

Polymer  Mechanics,  1965. 

50  Hill  R  ,  "A  Theory  of  the  Yielding  and  Plastic 

Hdtals-.  Proc.  Royal  Society,  Senes  A.  1948. 

Cl  r«io  R  V  and  Pioes.  R.  B.,  "Filamentary  Composite  Laminates 
^J«Ld  to  BiaxUI  sireU  Fields".  Tech.  Rep..  AFF0L-TR-73-n5. 
Air  Force  Flight  Dynamics  Lab.,  Wright-Patterson  Air  Force  Base, 

Ohio. 


78 


52  Bert,  C.  W.,  "Static  Testing  Techniques  for  Filament-Wound 
Coovosite  Materials'',  Composites,  January,  1974. 

T  M  and  Liber,  T.,  "Lamination  Residual  Stresses  In 
“•  NASA  CA-n5085.  June  1976. 

u  Petit  P.  H.,  "A  Simplified  fiethod  for  Determining  the  In-plane 
Shear *Stress-Stra in  Response  of  Unidirectional  Composites  , 

A.S.T.M.  STP  460,  1969. 

cc  dmah  r  U  "A  Simple  Procedure  for  Experimental  Determination 
“•  JinSUyoVnanJear  Nodylus  of  Unldiroctlonal  Coa,posuas". 

J.  Composite  Materials,  Oct.  1972. 

M  <;iiK  D  F  “In-plane  Shear  Stress-Strain  Response  of  Uni - 
directional  Composite  Materials",  J.  Composite  Materials,  Jan. 

1973. 

-  .  u  c  voftw  Y  T  "An  Investigation  of  the  Failure 

aIid"F?2ctIlre  Behavior  of ’oraphite/Epoxy  Laminates",  V.P.I.&  S.U., 
Report.  VPI-E-75-23,  September,  1975. 

ee  o.{.,cnn  H  F  and  Yeow,  Y.  T.,  "An  Experimental  Study  of  the 

Fracture  Behavior  of  Laminated  977 

CoBposite  Materials:  Testing  and  Design,  ASTM  STP  617,  1977. 

M  Yen-  Y  T  ,  and  Brinson,  H.  F.,  "The  Fracture  Behavior  of 

Gwhite/Epoxy  Laminates",  SESA  Spring  Meeting,  Dallas,  Texas, 

MV  15-20,  1977. 

60  Yeow  Y  T.,  and  Brinson,  H.  F.,  "A  Correlative  Study  Between 

gS'ii2sri--i.i-s5  ssr.r.:is:  ss:"-.. 

Vol.  15.  1967.  pp.  279-297. 

62  Kelly  A  ,  "Interface  Effects  and  the  Work  of  Fracture  of  a  Fibrous 
^oliteK  P?oc.  Roy.  Soc.,  A  319.  1970.  pp.  95-116. 

63.  Piggott.  M.  R..  "Theoretical  Estimation  of  Fracture  Toughness  of 
Fibrous  Composites",  J.  Mat.  Sci.,  Vol.  5,  1970,  pp. 

64.  Harrison.  N.  L.,  "The  Toughness  0^  « 

Co^Msite",  Fibre  Sci.  and  Tech.,  Vol.  4,  1971.  pp.  lOi  i'‘». 

65  Haddoups,  M.  E.,  Eisenmann,  J.  R., 

scooic  Fracture  Mechanics  of  Advanced  Composite  Materials  ,  . 

Coov.  Mat.,  Vol.  5,  1971,  pp.  446-454. 


p.  67. 

j  ••  r  D  **Pr*Artiur6  Anfllvsis  0^  Unidi rcc^iOHil 
67.  Sih,  G.  C..  and  .nan  pp!  230-2«. 

Composites  p  J.  Comp.  Mat.,  voi. 

pp.  343-359. 

..  /.  ru—  c  P  •  Huana  S.  L.;  and  McQuillen,  E.  N., 

69.  Sih,  G.  C.;  Chen,  E.  P.,  u  g,  ^  Filaner.t-Reinforced 

■Katerial  ^  IGh!  Sp.  167-126. 

Composites  ,  J.  Comp.  Hat.,  voi.  y,  hf 

Hat.,  Vol.  8,  1974,  pp.  253-265. 

a  D  1  anH  Uhitnev  J  M.,  "Uniaxial  Failure  of  Composite 

"■  Sinra;iaVnUl‘nrn,“  ;S  ^C^onceni^,tiona^ 

inertcan  Society  for  Testing  and  Materials.  197S.  pp. 

72  Cruse.  T.  A..  ’Tensile  Strength  of  Notched  Comosltes’,  J.  Coo. 
£?!!  Voi.  7.  1973.  pp.  218-229. 

”•  ?t"rl?gJh"in^6;tiSt;o;jJSiJe:4oxy^ 

American  Society  for  Testing  and  Materials,  iu/p, 

11,  1975,  pp.  315-328. 

and  Materials,  1973,  pp.  65-97. 

,  _  o  u  Fiiiicami  S  V  .  and  McLaughlin,  P.  V.,  Jr., 

”•  !?”;;4'aM‘F5tl9:e  M^chiniims  In  !!sME 

iwlastic  Behavio?  of  Cooosite  Materials.  AMD  -  Vol.  13.  ASMt. 

1975. 

77.  Kanninen,  M.  F.,  J^JJ^^^J^'^AnalysirAode/for  Crac^GJowth  in  a 
Development  of  a  ASTH-STP  617,  American 

Fiber  Reinforced  Composite  53-69 

Society  for  Testing  and  Materials,  1977,  pp.  53  69. 

”•  iSsSJoMUio;  6;’SjT?ni’a£?“|eoTG^^^^^  8oa«. 

Science  Publishers,  Inc.,  New  York,  1968,  p.  is/- 


80 


Vol.  13,  ASME,  1975.  pp.  73-96. 

•n  o™  u  F  and  Srawley,  J.  E.,  Plane  Strain  Crack  Toughness 

^  srKo  ‘ 

Force  Materials  Laboratory,  1974. 

H  1  Jr  "Mode  I  Stress  Intensity  Factors  for  Sym- 
82#  Konisri^  Mr«»  u  actm  ctp  Anipricsn 

Materials  Laboratory,  1973. 

M.  c™se,  T.  A.,  -user's  Manual  f>r  TAOJCRX:  A  2-D  Fracture  Mechanics 
Deck",  Pratt  and  Whitney  Aircraft,  1975. 

85  Oagoar.  S..  and  Broutman,  L.  J..  "Crack  Growth  Resistance  of 
Ran^  Fiber  Composites",  Journal  of  Composites,  Vol .  9,  July 

1975. 

86  Dally  J  W..  and  Alfirevich.  I.,  "Application  of  Birefringent 
Ratings  to  Glass-Fiber  Reinforced  Plastics".  Experimental 
Mechanics,  March  1969. 


APPENDIX  A 

STRESS-STRAIN  BEHAVIOR  AT  ROOM  TEMPERATURE 

The  purpose  of  this  section  is  to  present  an  analytical  tech¬ 
nique  for  deterwining  the  non-linear  stress-strain  response  to  failure 
of  syninetric  general  laminates  subjected  to  in-plane  loads  using  the 
lamina  non-linear  stress-strain  response  to  failure  properties.  In 
this  technique,  the  constitutive  relations  and  concepts  of  laminated 
plate  theory  [43.44]  are  employed  in  an  incremental  manner  to  predict 
the  behavior  of  general  laminates  under  monotonic  loading.  As  such, 
lamina-wise  stress  calculations  can  be  performed  and  sub-critical 
failures  can  be  predicted.  For  the  sake  of  brevity,  only  pertinent 
discussion  and  results  will  be  presented  here  as  more  extensive  docu¬ 
mentation  can  be  obtained  in  references  45  to  47.  In  addition,  an 
excellent  survey  and  discussion  on  the  types  of  analytical  techniques 
proposed  to  date  has  been  given  by  Rowlands  [48]. 

Method  of  Analysis 

A  non-linear  stress-strain  curve  can  be- approximated  quite  well 
by  a  series  of  linear  increments  provided  that  sufficiently  small 
increments  are  taken.  Assuming  that  an  increment  of  stress  is  pro¬ 
portional  to  an  increment  of  strain,  the  constitutive  relations  for  a 
lamina  under  generalized  plane  stress  for  a  given  increment  can  be 
expressed  as, 

**®i  *  ^ij^®i^  *  1»2.6) 
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that  Is, 


dOi^ 

^11  ^12  ^16 

dc,' 

da2 

►  3E 

C21  C22  ^26 

* 

dcj 

-^61  ^62  ^66- 

.'<”6 

or  simply, 

(do)  »  [C]  (de) 

where  do^.  de^  are  the  stress  and  strain  increments  respectively,  and 
the  stiffness.  C^y  is  a  function  of  Oy  In  the  case  of  a  general 
laminate,  the  incremental  constitutive  relations  for  the  lamina 
are  written  as 

{do*^}  -  CC*']  {dc'^}  (A-2) 

The  implication  of  equation  (A-2)  is  that  the  applied  stress  and 
material  axes  coincide.  Generally,  this  is  not  the  case.  However,  the 
two  coordinate  systems  can  be  related  by  a  transformation  matrix  [1,2], 
that  is, 

{do*'}  -  [T*^]  {d^} 

and, 

Idc'^}  »  [T*']  {di^>  (A-4) 

where  d^  and  dl*  are  the  stress  and  strain  increments  in  the 
lamina  with  respect  to  the  load  or  global  axes.  Substituting  equations 
(A-3)  and  (A-4)  into  (A-2)  yields, 

(d^)  ■  [C^]  (dc*^) 


(A-5) 
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where. 


[C*]  •  [t'']"’  [c'']  [t'*] 


tesinliq  that  the  stresses  are  uniform  through  the  thickness  of  the 
lamina,  the  resulunt  force  increment.  IdM).  with  respect  to  the  global 


coordinates  is  given  by, 


(dN)  «  z  t.Cdo  } 

k-1 


(A-6) 


where  L  is  the  thickness  of  the  k*"  lamina  and  L  is  the  total  n»d.er 
of  laminae  in  the  laminate.  Combining  equations  (A-5)  and  (A-b).  the 


r6su1t3nt  forcB  incrBroBnt  bccofUBS 


{dN}  -  z 

I _ ^  ^ 


(A-7) 


By  further  assuming 
laminae,  that  is. 


that  the  incremental  strains  are  the  same  in  all 


(dl*'}  -  {de*} 


(A-8) 


and  equation  (A-7)  reduces  to. 


where 


(dN)  «  [AlCde") 


k»L  1, 
[A]  -  r  tJJr} 
k«l  ^ 


(A-9) 


Equation  (A-9)  can  also  be  expressed  as, 

(dc*>  -  [A]"''(dN} 


(A-10) 


In  equation  (A-10),  [A]"’  represents  the  average  laminate  compliance 
during  the  (n+1)^*’  load  increment.  However,  these  properties  are  not 


84 


known  when  the  load  increment  is  applied.  Consequently,  the 

princlp.1  elastic  properties  at  the  end  of  the  n«'  load  increnent  are 
used  to  compute  the  laminate  compliance  for  the  (ntl)“  load  increment. 
Sequentially,  knowing  the  (ntl)^''  load  increment  and  the  laminate 
compliance  at  the  end  of  the  n‘''  load  increment,  the  laminate  strain 
increment  idc")  and  each  lamina  strain  increment  {dc*")  and  stress  incre- 
nent  (d^l  can  be  calculated  by  using  equations  (A-10).  (A-4)  and  (A-2), 
respectively.  These  laminae  stress  and  strain  increments  are  added  to 
the  stresses  and  strains  at  the  end  of  the  n"  load  increment  to  ootain 
the  current  stresses  and  strains  in  each  of  the  laminae. 

The  incremental  loading  procedure  described  above  is  continued 
until  a  lamina  reaches  a  state  of  stress  when  it  can  no  longer  sustain 
additional  load.  Analytically,  this  state  of  stress  is  detem-ined  by 
a  failure  criterion.  Numerous  failure  criterion  are  available  [27] 
for  predicting  the  strength  of  composite  materials  and  these  criteria 
can  be  classified  as  either  having  independent  or  dependent  failure 
«Kles.  The  maximum  stress  or  strain  criteria  are  the  ones  with 
Independent  failure  modes.  The  criteria  proposed  by  Ashkenazi  [49], 
Hill  [50]  and  Puppo  and  Evensen  [34]  are  examples  of  failure  criteria 
with  dependent  failure  modes.  Between  the  two  classes,  experimental 
results  [51]  tend  to  confim  that  the  criteria  with  dependent  failure 
■odes  are  more  appropriate  for  composite  materials. 

The  assumptions  made  in  Ashkenazi’s  criterion  are  that  the 
material  is  a  macroscopic  continuum,  the  strength  properties  are 
fourth  order  tensors  and  environmental  effects  are  neglected.  His 
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criurton  for  a  pUno  orthotropic  material  onPer  oniaxial  load  it. 

1_ .  CosS  .  (J.  '  11  cot'e  Sio^e  t  51^  (H-'D 

Ojj  X  1X45  * 

^re  0,  it  the  applied  axial  ttrett.  e  it  the  angle  between  0,  and 
the  direction  of  the  applied  load  and  X,  Y  and  X„  are  the  tentile 
strengtht  along,  trantverte  and  at  45”  to  the  fibert,  respectively. 
Mile  eguatioo  (A-ll)  is  written  for  uniaxial  loading  only,  other 

states  or  stress  were  considered  by  Ashkenazi. 

using  assomptlont  similar  to  Ashkenazi.  Hill  developed  an  aniso¬ 
tropic  yield  criterion  which  can  be  modified  into  a  failure  criterion 
for  an  orthotropic  material  under  a  generalized  state  of  plane  stress. 

It  can  be  expressed  as« 


X 

2 

®2 

^  °1®2  ^ 

^12^ 

+ 

Y 

XY 

h  4 

k  * 

M«re  the  stresses  o,.  o^  and  are  along  the  material  axes.  S  is  the 
intral»nina  shear  strength  and  X  and  Y  are  as  defined  for  equation 

(A-11). 

Finally,  a  more  general  failure  criterion  proposed  by  Puppo  and 
Evensen  can  be  expressed  as 


(A-13) 
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^ere  ,  •  ^  H  >"  <"«raction  factor  aod  the  other  quantities  are  as 
previousl/defined.  Further,  the  authors  indicated  that  their  theory 
could  he  adaoted  to  a  -ide^range  of  naterials  b,  writing  the  inter¬ 
action  factor  as  r  •  (|f^"  the  exponent,  n.  is  a  material 

paraneter. 

so  far,  the  discussion  has  been  related  to  the  detemination  of 
the  states  of  stress  In  a  lamina  and  the  subsequent  failure  resulting 
therefrom.  Tl»  most  difficult  and  fundamental  question  to  answer  in 
such  inuestigations  Is  how  to  represent  failed  laminae  mathematically 
uhen  s«  lairtnae  in  the  laminate  are  still  intact  or  hare  not  reached 
their  critical  stress  state.  An  easy  analytical  technique  is  to 
mathematically  remre  failed  laminae  from  the  laminate  by  zereing  the 
thicknesses  arwl  stiffnesses  of  the  failed  laminae.  Howerer.  this  pro¬ 
cedure  leads  to  conservative  or  ultra-conservative  strength  predictions 
depending  on  the  laminate  lay-up.  An  alternative  is  to  uncouple  the 
constitutive  equation  of  the  failed  laminae  from  the  laminate  stiff¬ 
ness  or  coniliance  calculations.  This  can  be  done  by  assigning  zero 
stiffnesses  to  the  three  off-diagonal  terms  of  the  matrix  in 
equatior  (A-5)  and  the  appropriate  diagonal  term  when  its  ultimate 

Strain  is  reached  or  exceeded. 

Results  and  Discussion 

Using  the  incremental  procedure  previously  explained,  the  total 
stress-strain  response  of  general  laminates  subjected  to  uniaxial  load- 
ing  are  generated.  The  results  obtained  are  shown  in  Figures  A-1  to 
A-3.  figure  A-1  shows  the  analytical  and  experimental  results  for  the 


Figure  A-1. 
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[0“/i30“/0®325  laminate.  Excellent  correlation  is  seen  up  to  about 
50  ksi  (344.8  MPa).  At  this  stress  level,  failure  of  the  ±30®  laminae 
is  predicted  and  is  indicated  by  the  steeper  slope  of  the  experimental 
stress-strain  curve.  After  this  sub-critical  failure  stress  level,  the 
analysis  predicted  higher  stresses  than  experiment.  This  is  probably 
due  to  the  omission  of  interlamina  stresses.  Reasonable  correlation 
betueen  analysis  and  experiment  is  observed  in  the  case  of  the 
[45*/15*/75"/45“]25  laminate.  Figure  A-2.  Excellent  correlation  is 
seen  in  the  case  of  the  [90®/t60®/90®]25  laminate.  Figure  A-3. 

A  comparison  of  the  normalized  laminate  failure  strengths,  as 
predicted  by  the  three  failure  criteria  used  in  the  analysis,  are  as 
shown  in  Table  A-1.  Among  the  three  failure  criteria  used,  Ashkenazi’s 
appears  to  be  the  best.  It  should  be  noted  that  the  exponent  of  the 


Table  A-1.  Predicted  Failure  Strength,  Normalized  with  Respect  to 
the  Experimental  Strength. 


Laminate 

Orientation 

Ashkenazi 

Hill 

Puppo-Evensen 

[(r/i30*/0®]25 

0.971 

0.965 

0.968 

[45V15V75V45®]25 

1.002 

1.168 

1.054 

[90*/±60*/90»]2s 

0.827 

0.786 

0.864 

interaction  factor  in  equation  (A-IS)  was  taken  as  unity.  By  adjusting 
this  exponent,  for  the  laminates  investigated,  identical  analytical  and 
experimental  strengths  can  be  obtained.  Unfortunately,  the  exponent 
MS  not  constant  even  though  the  specimens  were  obtained  from  a  single 


large  panel.  Nevertheless,  the  Puppo-Evensen  criteria  appear  to  be 
the  most  promising  as  it  can  be  used  for  multi-axial  loading,  allows 
for  material  variations  through  the  interaction  factor  and  can  be  used 

without  a  lamina-wise  analysis  as  shown  here. 

On  the  whole,  good  stress-strain  response  was  obtained  from  the 

non-linear  Incremental  analysis  proposed  here. 

Application  of  Analysis 

Besides  the  toUl  stress-strain  prediction  capability,  the  incre- 
iienul  procedure  described  here  can  be  used  to  validate  test  methods 
such  as  intralamina  shear  response  determination.  Presently,  a  variety 
of  test  methods  has  been  proposed  [52]  to  obtain  this  in-plane  shear 
stress-strain  response.  Because  among  the  four  principal  material 
properties  reguired  in  the  two-dimensional  analysis,  the  most  difficult 
esperimenul  quantity  to  obuin  is  the  intralamina  shear  stress-strain 
response.  This  difficulty  is  best  clarified  by  Figure  A-4.  where 
variations  in  the  response  of  five  different  test  methods  are  observed. 
In  this  figure,  the  10-deg  off-axis  test  method  of  Daniel  [53]  and 
Chamis  et  al  [35]  and  the  calculated  [15]3,  results  appear  to  correlate 
In  the  case  of  the  tensile  test,  the  methods  of  Petit  [54]  and 

Rosen  [55]  tend  to  agree.  The  results  obtained  from  the  [OVSO’],^ 
syrotric  rail-shear  method  of  Sims  [51]  indicates  the  most  ductility. 
Thus,  depending  on  the  test  method  used,  the  intralamina  shear  stress- 
strain  response  may  vary  from  a  brittle  to  ductile  nature. 

To  validate  the  various  test  methods,  an  easy  method  is  to  com¬ 
pare  the  actual  initial  axial  modulus  of  the  off-axis  tensile  tests 


ermine 


with  the  orthotropic  transfomation  equation  which  can  be  expressed 


as, 


2v 


12 

n 


2  2.  .  Sin^e 

Cos  e  Sin  e  +  -p — 

42 


(A-U) 


w«l  in  equation  (A-14)  are  as  previously  defined.  By 
substituting  the  initial  value  fron  tbe  off-axis  calculattd 

,^\a  and  the  initial  G,^  value  by  Rosen's  [t45»]j  method  into  equa- 
tim  (A-14).  an  analytical  and  experimental  conparison  can  be  made  as 
Figure  A-5.  In  this  figure,  excellent  correlation  is  seen 
b.t>«en  experiment  and  analysis  when  the  [15'']85  results  are  used.  It 
ut  be  noted  that  the  reason  for  the  choice  of  these  two  methods 
amng  the  five  methods  is  that  the  responses  of  the  two  chosen  methods 

appear  to  be  the  average  upper  and  lower  bounds. 

To  further  validate  the  effect  of  the  two  lamina  shear  responses 
on  laainate  behavior,  the  total  response  curves  obtained  from  the  two 
Kthods  (ClS-las  calculated  and  Rosen’s  [t45»]^  method)  are  used  in  the 
increMital  analysis,  previously  described,  with  the  other  three 
parameters  (E,,.  and  v,2) 

behavior  of  general  laminates.  The  analytical  results  obtained  from 
the  two  shear  results  are  compared  to  the  experimental  tensile  test 
resulU  as  shown  in  Figures  A-6  to  A-8.  In  the  case  of  the  [0V90*]4s. 
[0-/±45*/0*]2s.  [90“/±60“/90'*]23  laminates.  Figures  A-b  to  A-7. 

little  difference  between  the  analytical  results  were  observed. 
Apparently  in  these  cases.  Intralaminar  shear  response  does  not  play 
an  iiportant  role  as  in  the  case  of  the  [45V0V90V-45‘*]2s. 


Tensile  Modulus 


(^S5|  -01  *)  smnaow  3iisN3i 


FIBER  ORIENTATION  (degrees) 

Variation  of  the  Tensile  Moduli  with  Fiber  Orientation  for  Uni 
directional  T300/934  Graphite/Epoxy  Laminates. 


Experimental  Results 

[15'']8s  Off-Axis  Shear  Behavior 

[i45“]^5  Shear  Behavior 
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Experimental 

[15"]qj  Off-Axis  Shear  Behavior 
[t45“]^.  Shear  Behavior 
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.nd  [60-]j,  laminates  shown  In  Fi,ures  A-7  to  A-S.  -hero  ?o«d  correla¬ 
tions  are  obtained  from  the  [U-lgs  calculated  results. 

The  main  conclusion  drawn  from  these  results  is  that  the  off-axis 
test  is  the  better  method  for  the  detennination  of  the  shear  response 
of  unidirectional  laminates.  On  the  other  hand,  the  [=45=]^  tensile 
and  [0“/90'*]j  rail  shear  methods  seem  to  include  interlamina  effec 
Thus,  it  applars  that  the  apprppriate  test  method  would  depend  primaril, 
upon  what  is  sousht.  For  example,  if  the  effect  of  a  variable,  sa, 
environment  (temperature  and  humidity),  on  the  shear  properties  wsthin 
a  lamina  are  to  be  detemined,  the  off-axis  test  would  be  appropriate. 
If.  on  the  other  hand,  the  effect  of  the  same  variables  on  laminatt 
shear  response,  including  Intralamina  and  interlamina  effects,  the 
[145*],  tension  or  the  [OVOO'lj  rail  shear  test  would  be  the 

appropriate  test. 


APPENDIX  B 


A  SYNOPSIS  OF  FRACTURE  BEHAVIOR  OF  NOTCHED  COMPOSITES 
AT  ROOM  TEMPERATURE 

In  this  appendix,  a  brief  perspective  of  some  of  the  proposed 
analytical  models  for  the  fracture  behavior  of  notched  composites  is 
presented.  This  is  followed  by  the  analytical  methodology  used  for  the 
nacroscopic  fracture  behavior  of  composites,  experimental  considerations, 
and  an  evaluation  of  the  available  fracture  models  with  available 

experiaental  data.  The  application  of  these  concepts  are  used  to 

discuss  possible  future  time-temperature  investigations  on  notched 
conposites.  Only  pertinent  discussion  and  results  are  presented  here 
as  more  extensive  documentation  can  be  obtained  in  references  57  to  60. 

Fracture  Model s 

The  fracture  behavior  of  composite  materials  has  been  studied 
using  mainly  three  aspects:  Linear  elastic  fracture  mechanics  (LEFM) 
on  a  microscopic  scale  with  the  material  being  heterogeneous  and 
anisotropic;  LEFM  on  a  macroscopic  scale  with  the  material  being 
hoiogeneous  and  anisotropic;  and  the  material  modelling  approach  using 
approximate  models  to  represent  the  major  effects  of  heterogeneity  and 

anisotropy. 

Chronologically,  the  fracture  behavior  of  notched  composites 
was  analyzed  from  the  micromechanical  viewpoint  [61-64]  first  and  later 
was  exwined  from  the  macromechanical  viewpoint  [65-74].  The  primary 
limitations  of  the  micro  and  macro  approaches  are:  The  former,  due  to 
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ns  siwiicttj.  is  not  feasible  -hen  n«Ui-axial  loadine  is  considered 
and  in  the  latter,  the  non-self-similar  crack  propagation  (which  exists 
for  certain  laminate  lay-ups)  is  not  considered.  Since  the  mathe- 
Mtlcal  theory  of  LEFH  assumes  self-similar  crack  propagation,  the 
latter  appreach  will  not  be  applicable  when  this  assmrption  is  .iolattd. 
To  o.erco«  these  limitations,  the  “Materials  Science-  model  [75,76] 
and  the  “Local  Heterogeneous  Region"  (LHR)  model  [77]  were  proposed. 

The  adyantages  gained  in  these  material  modelling  methods  are  that  the 


structural  response  and  various  failure  modes  and  mechanisms  can  be 
incorporated  and  predicted  regardless  of  the  boundary  loadings.  Conse 
guently.  the  n»*er  of  parameters  reguired  by  the  models  increases  as 

its  analytical  capabilities  increase. 


Analytical  Methodology 

The«de  I  critical  stress  intensity  factor.  Me.  for  an  infinite 
isotropic  plate  containing  a  crack  Is  expressed  as 


(B-1) 


and  the  nodified  Kj^ 
expressed  as 


using  an  Irwin  type  plasticity  correction  factor  is 


o^/xla  +  a^) 


(B-2) 


where  is  the  critical  remote  stress  normal  to  the  crack,  2a  (Figure 
B-1)  is'the  toUl  crack  length,  and  a,  is  the  characteristic  length  or 
inherent  flaw  length  as  denoted  by  Waddoups  et  al.  [65].  These  authors 
assumed  that  an  unnotched  laminate  is  equivalent  to  a  notched  laminate 
with  a  notch  length  of  2a^.  That  is.  for  the  unnotched  laminate 
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0  0 


(8-3) 


wher«  a  is  the  ultimate  strength  of  the  unnotched  laminate.  By  per¬ 
forming  txo  tests,  one  on  an  unnotched  laminate  and  the  other  on  a 
notched  laminate,  the  two  parameters  a^  and  are  determined.  Using 
equations  (B-2)  and  (B-3),  the  critical  remote  stress,  c^,  for  the 
same  laminate  containing  another  flaw  size,  2b,  is  related  to  and 

a^  as  [65] 

!o.  Fin  (B-4) 


Using  a  similar  procedure  as  Waddoups  et  al .  but  with  a  different 
rationale,  Whitney  and  Nuismer  [70,71]  use  the  exact  anisotropic 
center  cracked  solution  of  Lekhnitskii  [73]  to  derive  two  failure 
critena  for  center  notched  specimens.  The  normal  stress  distribution 
for  an  infinite  plate.  Figure  B-1,  containing  a  crack  and  subjected  to 

uniaxial  tension  o  is 

/  1  -  (8-5) 

"  42  /Tta(x2  -  a^ 

where  Kj  from  LEFM  is  expressed  as  Kj  =  o/^a.  The  first  failure 
criterion  of  Whitney  and  Nuismer,  known  as  the  point  stress  criterion. 
States  that  failure  of  a  notched  laminate  will  occur  when  the  local 
stress  0  (a^.o)  at  a  certain  distance,  a^,  from  the  implanted  notch  tip 
reaches  the  strength  of  the  unnotched  laminate,  o^.  Using  equation 
(B-5)  and  the  definition  of  the  point  stress  criterion,  results  in 
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(B-6) 


where 


a  ♦  a. 


Ttw  second  failure  criterion,  known  as  the  average  stress  criterion, 
suted  that  failure  will  occur  when  the  average  stress  o^(d^.o).  over 
a  distance  d^  from  the  notch  tip  reaches  the  strength  of  the  unnotched 
laminate,  oj  Substituting  equation  (B-5)  into  the  definition  of  the 
average  stress  criterion  gives 


“c .  /1Z5 

vrri: 


(B-7) 


where 


•1 


TT 


The  paraneters  o^.  a,,  and  d^  are  found  experimentally  using  the  same 
two  required  tests  as  in  the  Waddoups  et  al.'s  method.  As  discussed  by 
Tsai  and  Hahn  [79],  the  characteristic  lengths  d^^  and  a^  are  related  by 

a  factor  of  two,  that  is,  d^  * 

As  shown  in  reference  [70],  the  fracture  toughness  according  to 
the  point  stress  criterion  is  given  by 

liq  ■  V»«(’  - 

whorus  for  tho  average  stress  criterion  the  fracture  toughness  is 


Kg  ■  Og/naO  - 


(B-9) 
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Th.  analysis  discussed  so  far  is  for  an  infinite  plate  while  for 
econonic  reasons  and  experimental  convenience  finite  width  plates  are 
used.  In  order  to  correlate  the  two  types  of  results,  the  experimenul 
dau  is  corrected  to  the  analytical  results  by  applying  appropriate 
finite  width  correction  factors.  The  usual  procedure  is  to  detemine 
the  infinite  plate  notched  strength  from  the  equation 

o,  .  F(a/ll)o/ 

,t«re  Oj  is  the  infinite  plate  notched  strength.  0/  is  the  notched 
strength  of  a  finite  width  specimen,  and  F(a/Vi)  is  the  finite  width 
correction  factor.  The  value  of  F(a/ll)  is  determined  from  the  eguation 

for  isotropic  materials  [80]. 


F(a/W)  »  /(W/«a)  taniwa/w) 


(B-11) 


where  W  is  the  width  of  the  finite  plate. 

Tl*  analytical  methodology  mentioned  to  this  point  depends  very 

strongly  on  the  isotropic  treatment  of  the  analysis  as  seen  in  the  cor- 
relation  of  the  finite  width  experimental  data  with  the  infinite  plate 
auilytical  solution,  or  vice  versa.  However,  the  difference  between  the 
correction  factor  for  isotropic  and  anisotropic  materials  is  of  the 
order  of  10!,  as  shown  by  Cruse  and  Osias  tSl]  and  Konish  [823.  Thus, 
as  a  matter  of  convenience,  the  isotropic  correction  factor  is  used. 

To  avoid  the  assumption  of  an  isotropic  factor.  Snyder  and  Cruse 
[,4.833  developed  a  boundary-integral  equation  (BIE)  method  to  calculate 
finite  width  effects  and  stress  Intensity  factors  for  modes  I  and  11. 
for  various  geometries  and  in-plane  loadings.  In  this  approach  the 
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■aterial  is  rodelled  as  a  homogeneous,  linearly  elastic,  anisotropic 
continue  using  the  assumptions  of  lamination  theory.  The  BIE  method  is 
based  upon  an  appropriate  fundamental  solution  to  the  field  equations 
of  elasticity.  The  fundamental  solution  is  used  to  derive  a  set  of 
integral  equations  relating  boundary  values  of  displacement  and  traction 
corresponding  to  a  well-posed  boundary  value  problem.  The  boundary 
integral  equations  are  solved  using  approximate  techniques,  e.g.. 

■Ddeling  the  boundary  by  constant-displacement,  constant- tract! on 

segments. 

The  boundary  integral  equations  relate  boundary  values  of  traction 
and  displacenent  on  that  portion  of  the  boundary  excluding  the  crack 
boundary,  thus  avoiding  the  problem  of  modeling  the  crack  boundary. 

An  additional  advantage  of  the  method  is  that  it  provides  a  direct  means 
for  obtaining  both  of  the  stress  intensity  factors  Kj  ird  Kjj.  i.e.. 
equations  are  derived  for  obtaining  Kj  and  Kjj  directly  from  the 
boundary  solution.  Thus,  the  boundary  solution  may  be  computed  using 
other  n«erical  methods,  such  as  the  finite  element  method. 

As  developed  by  Snyder  and  Cruse,  the  basic  equation  of  the  BIE 

method  may  be  written  as 

Uj(P)/2  +  I  u^{Q)  Tj^(P.Q)ds(Q)  »  Uj.{P.Q)ds(Q)  (B-12) 

In  equation  (B-12).  u^(Q)  and  t^(Q)  refer  to.  respectively,  displacements 
and  tractions  at  boundary  points  Q,  not  including  the  crack  boundary. 

The  BIE  method  uses  special  kernels,  Tj^(P,Q)  and  Uj-(P,0),  for 
Tj^(P.Q)  and  Uj.(P.Q)  in  equation  (B-12)  to  mathematically  model  a 
traction-free  straight  crack  of  length  2a.  Since  T^.  and  t.  are 
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identically  zero  on  the  crack  surface,  the  crack  need  not  be  modelled. 

The  expression  for  determining  mode  I  and  mode  II  stress  intensity 
factors  is  an  integral  of  boundary  displacements  and  boundary  tractions 
on  an  arbitrary  path  in  the  body,  and  is  given  by 

Kr  ft  *  f  R!*"(Q)u.(Q)ds(Q)  +  f  Lj‘“(Q)t.(Q)ds(Q)  (B-13) 

I, II  1  '  Js 

The  kernels  and  are  complex  variable  functions  which  include 
crack  size  and  location,  as  well  as  material  properties.  As  discussed 
in  referwce  [84],  while  the  derivation  of  equation  (B-12)  leads  to  the 
definition  of  equation  (B-13),  equation  (B-13)  may  use  boundary  data 
taken  from  any  numerical  model  of  the  cracked  geometry.  The  only  input 
infonnation  needed  in  the  technique  is  geometry,  material  properties, 
and  grid  sizes  and  spacings. 

Experimental  Considerations 

In  order  to  determine  the  versatility  of  the  three  two-parameter 
models  and  the  one  one- parameter  model  described  in  the  previous  section 
a  series  of  center  notched  specimens  containing  narrow  slits  with 
variations  in  aspect,  ratio  (2a/w.  where  w  is  the  total  width  of  the 
specimen)  and/or  anisotropy  were  tested.  In  the  second  series  of  tests, 
the  aspect  ratio  variations  were  the  same  as  the  first,  but  the  flaw 
geometry  was  changed  to  that  of  a  hole  or  square  slot  rather  than  a 
slit.  The  objective  of  the  second  series  was  to  ascertain  the  types  of 
fracture,  if  more  than  one  existed  and  the  fracture  strength  variation 
when  the  height  of  the  slit  was  increased.  In  these  series  of  tests, 
experimental  methodologies  such  as  crack  opening  displacement  (COO) 
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gages.  t>«  types  of  photoelastic  coatings  and  photographic  recording 
of  crack  propagation  were  employed.  In  addition,  in-situ  and  post 
.icroscopic  observation  of  the  crack  growth  and  similar  observations 
of  the  free  surfaces  of  the  specimens  were  made. 

In  the  third  series  of  tests,  the  implanted  flaw  geometries  and/ 
or  location,  material  anisotropy  and  aspect  ratios  were  again  varied. 

The  notch  geanetries  and/or  locations  chosen  were  the  single-edged 
notch  (SOI),  double-edged  notch  (DEN)  and  centrally  located  circular 
holes.  The  reasons  for  this  series  of  tests  were  to  compare  the 
fracture  toughnesses  of  the  SEN  and  DEN  specimens  using  the  BIE  method 
and  the  eiperimental  fracture  strengths  of  the  specimens  containing 
holes  with  the  predictions  of  the  three  two-parameter  models. 

Results  and  Discussion 
COO  Measurements 

Typical  load-COD  traces  for  all  the  laminates  in  the  first  series 
of  tests  .«  as  show,  in  Ftgute  B-2.  The  wasureinents  -ere  taken  fro. 
the  relatite  lotioo  of  the  center  of  the  opposite  faces  of  the  notch 
In  the  direction  of  the  remote  tensile  load.  Similar  traces  were  also 
obtained  from  the  second  series' of  tests  (slots).  Generally,  as  indi¬ 
cated  in  Figure  B-2,  the  COVsdSVO'ljj  laminates  exhibited  nearly 
linear  response  while  the  and  laminates  were  highly 

non-linear.  In  most  cases,  the  displacements  of  the  two  latter  lami- 

nates  exceeded  the  range  of  the  COO  gage. 

Gaggar  and  Broutman  [86]  attributed  such  non-linear  behavior 

and  discontinuities  in  load-COD  measurements  to  defonnations  and 


CRACK  OPENING  DISPLACEMENT  (x  lO"  cm) 


displacement  curves 
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extensions  of  the  notch  tips.  Assuming  self-similar  cracit  propagation, 
those  authors  suggested  a  compliance  matching  procedure,  which  is 
illustrated  in  Figure  B-2,  to  quantify  the  amount  of  actual  or  apparent 
crack  extension  at  any  arbitrary  load  level.  Using  this  procedure, 
the  compliances  of  the  slots  were  always  higher  than  their  slit 
counterparts.  The  more  flexible  nature  of  the  slots  was  expected  as 
there  was  less  constraint  than  in  the  case  of  the  slits. 

Fracture  Strengths 

The  variation  of  the  average  value  of  the  critical  gross  stresses 
with  aspect  ratios  for  the  first  and  second  series  of  tests  (i.e.. 
slits  and  slots)  are  as  shown  in  Figure  B-3.  Although  the  trends  of 
the  critical  gross  stresses  were  similar,  the  slots  always  fractured 
at  a  higher  stress  level  than  the  slits  irrespective  of  the  aspect 
ratio  or  laminate  orientation.  This  was  probably  due  to  the  more 
flexible  nature  of  the  slots  as  determined  by  the  COD  measurements. 

In  the  case  of  the  laminate,  the  specimens  containing 

circular  holes  fractured  at  about  the  same  or  slightly  lower  stress 
levels  than  the  ones  containing  narrow  slits  for  all  aspect  ratios  as 
shown  in  Figure  B-3.  Examination  of  the  photoelastic  fringes  of 
Figures  B-4a  and  B-4c  tended  to  indicate  that  the  size  of  the  region 
of  the  intense  stress  gradient  was  larger  for  the  hole  than  for  the 
slit  at  about  the  same  remote  stress  level.  This  would  imply  that 
the  holes  were  more  critical  than  the  slits  for  the  [i45  3^^  laminate. 

It  should  be  noted  that  the  trends  of  the  variation  of  the 
with  aspect  ratio  were  the  same  for  all 


critical  reinote  stresses 


O.z  U.J  U.H 
2fl/w 

flour.  B-3.  noTOt.  F.nur..  Str.i..l  0»  1300/934  Gr.phit./Epoxy  U.mlMtox 


a)  13.10  ksl 
(90.29  HPa) 


b)  Fracture 


c)  13.22  ksl  d)  Fracture 

(91.16  MPa) 


e)  16.39  ksl  Fracture 

(113.00  MPa) 


B-4.  Isochromatics  of  [±45"]4s  Laminate  with  Perforated 
(a-d)  and  Continuous  (e.t)  Coatings  ..r./w  -  0.45). 


112 


cutouts  as  shown  in  Figure  B-3.  This  observation  tends  to 
subsuntiate  the  “Materials  Science"  model  for  fracture  which  assumes 
that  flaw  geometry  does  not  have  a  major  influence  on  the  laminate 
fracture  strength  and  the  strength  is  more  closely  related  to  the  size 
and  shape  of  the  damage  zone  surrounding  the  flaw  tip. 

Finally,  the  critical  remote  stresses  of  the  specimens  contain¬ 
ing  SER,  DEN  and  holes  with  aspect  ratio  are  shown  in  Figures  B-5, 

B.6  and  B-7.  respectively.  Except  for  the  laminate,  the  trends 

of  the  variation  of  the  fracture  strengths  with  aspect  ratios  in 
every  case  were  similar  to  the  ones  observed  for  isotropic  materials. 
The  deviations  of  the  laminate  from  these  trends  were  likely 

due  to  the  axial  crack  propagation  which  propagated  into  the  gripped 
sections  of  the  specimen  prior  to  transverse  failure. 

Birefringent  Coatings 

Typical  results  obtained  for  each  laminate  orientation  using  both 
continuous  and  perforated  coatings  are  shown  in  Figures  B-4c  and  B-4e. 
That  is.  the  fringe  density  was  less  in  the  case  of  the  continuous 
coatings  but  the  shapes  of  the  patterns  were  similar. 

Figure  B-8  shows  the  types  of  stress  fields  and  fracture  planes 
obUined  for  specimens  containing  center  slits.  In  the  case  of  the 
[0*]g  and  [0"/90“]^j  specimens,  although  both  laminates  could  be  con¬ 
sidered  to  be  orthotropic  and  the  intense  shear  stress  regions  were 
in  the  axial  direction  (see  Figure  B-8a  and  B-8b).  the  final  fracture 
planes  were  different.  Non-self-similar  crack  propagation  was 
observed  for  the  [O^as  laminates  and  self-similar  crack  propagation 
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MS  Observed  for  the  [OVSO"]^^  laminates.  The  intense  shear  stress  . 
regions  for  the  laminates  were  adjacent  to  the  notch  tips 

and  were  in  the  AS"  directions.  Generally,  the  pattern  was  more 
dense  in  the  +45‘>  direction  than  the  -45»  direction.  Since  the  outer 
laminae  were  in  the  +45°  direction,  this  observation  merely  indicates 
that  the  inner  -45°  laminae  were  constrained  by  the  outer  laminae 
(♦45°  direction)  to  produce  this  effect.  The  fringe  pattern  was 
similar  to  the  ones  observed  in  isotropic  materials  in  the  case  of  the 
[0Vi45V0°]2s  laminate;  i.e..  butterfly  wing  shaped  patterns. 

Figure  B-8d. 

No  attempts  were  made  to  quantify  the  birefringence  results  as 
the  mismatch  between  the  properties  of  the  coating  and  the  underlying 
mateHal  (Poisson's  ratio  and  moduli)  as  well  as  anisotropy  of  the 
laminate  made  proper  fringe  interpretation  difficult  [86].  Further, 
the  properties  of  the  material  were  not  constant  and  changed  with 
stress  level.  Similar  results  have  been  reported  by  others  [12]. 
Perhaps  more  importantly,  the  uncertainties  relative  to  the  nature  and 
strength  of  the  singular  stress  field  for  the  slits  and  slots  further 
compounded  our  fringe  interpretation  difficulties.  That  is.  the 
nature  of  the  singularity  of  the  laminate  was  likely  quite  different 
than  the  nature  of  the  singularity  in  either  the  continuous  or  perfo¬ 
rated  coating. 

In-Si tu  and  Post  Failure  Examinations 
Stable  crack  growth  was  often  observed  during  most  of  the  tests, 
In  all  cases  stable  crack  growth  occurred  by  matrix  fracture.  In  the 
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C.S.  of  the  CO"]8,  lo-inates,  steble  axial  crack  growth  alwa,t  occurred 
at  the  notch  tip  and  proceeded  to  the  gripped  sections  of  the  specinen. 
Catastrophic  transverse  failure  would  occur  only  upon  completion  of 
the  axial  crack  propagation. 

The  most  interesting  situations  of  stable  crack  growth  was  ob¬ 
served  in  the  DEN  [g0Vt60V90"]2s  tests.  A  typical  example  of  crack 
propagation  for  this  notch  geometry  is  shown  in  Figure  B-9.  Initially, 
««  self-similar  crack  would  originate  at  the  notch  tips  and  propaga¬ 
tion  would  continue  as  the  remote  load  was  increased.  At  higher  stress 
levels,  micro-cracking  was  observed  away  from  the  plane  of  the 
iwlanted  notch.  Finally,  separation  of  the  specimen  occurred  at  the 

region  where  the  cross-sectional  was  minimum. 

Post  failure  microscopic  examination  of  all  the  unnotched  and 
notched  specimens  of  the  first  and  second  test  series  revealed  a  large 
aater  of  micro-cracks  on  all  the  free  edges.  In  the  notched  cases, 
microscopic  examination  clearly  indicated  more  pronounced  micro-cracks 
on  the  free  surfaces  of  the  implanted  notch  as  opposed  to  the  specimen 

boundary  free  surfaces. 

A  qualitative  comparison  of  the  boundary  free  edges  for  unnotched 
and  notched  specimen  within  each  general  laminate  lay-up  showed  more 
pronounced  cracks  in  the  notched  case  than  in  the  unnotched  case.  This 
was  true  even  though  the  remote  stress  levels  of  the  unnotched  speci¬ 
mens  were  the  higher  of  the  two  situations.  This,  perhaps,  indicated 
that  the  notches  had  an  adverse  affect  on  the  stress  sute  at  the 
specimen  free  edge. 


low  the  notch  plane. 
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Characteristic  Lengths 

In  order  to  obtain  favorable  correlations  between  the  two- 
paraeter  fracture  models  (inherent  flaw,  point  stress  and  average 
stress)  and  experiments,  the  characteristic  length  or  damage  tone 
«Kt  be  ktwnm.  In  fact,  even  to  use  these  methods  at  all  requires 
that  this  dimension  must  be  established  for  a  particular  laminate. 

The  oHginators  of  the  inherent  flaw  and  point  stress  criteria  sug- 
,ested  a  characteristic  length  of  0.04  inch  (1.016  am)  for  all  aspect 
ratios  would  be  a  reasonable  assumption  while  for  the  average  stress 
criterion  they  suggested  a  characteristic  length  of  0.15  inch  (3.810  ™) 
They  further  suggested  that  this  quantity  could  be  assumed  to  be 
coRSUnt  for  a  type  of  composite;  i.e..  graphIte/epoxy.  glass/epoxy, 
etc.  To  validate  the  above  suggested  characteristic  lengths  for  our 
materials  (T300/934  graphite/epoxy),  an  attempt  was  made  to  calculate 
B«e  parameters  (Equations  (8-4),  (8-6)  and  (8-7))  from  the  experi- 

■ental  data  generated  in  this  investigation. 

As  Shown  in  Figure  8-10,  the  characteristic  lengths  determined 
froi  our  experirmnts  were  higher  than  the  ones  suggested  previously 
[66.72].  It  should  be  noted  that  non-self-similar  crack  growth  was 
Obsirved  in  the  cases  shown  in  Figure  8-10  and  it  is  probably  due  to 
this  fact  that  a  decreasing  characteristic  length  with  Increasing 
aspect  ratio  trend  was  observed.  The  proposed  two-parameter  fracture 
^is  do  not  allow  tor  such  a  trend!  In  all  fairness,  however,  the 
aawlels  were  originally  proposed  for  quasi-isotropic  laminates  and  the 
lailnate  types  shown  in  Figure  B-10  are  not  quasi-isotropic. 


2a  (im) 
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Quasi-self-similar  crack  growth  was  observed  for  the  nearly 
quasi-isotropic  laminates  and  the  trends  were  similar  to  the  ones  sug¬ 
gested  by  the  proposed  models  as  evident  by  examination  of  Figure  B-11. 
In  this  latter  figure,  the  total  crack  extension  prior  to  separation 
as  determined  by  the  load-COO  method  were  consistent  with  the 
characteristic  lengths  of  the  inherent  flaw  model  for  the  [0Vi45V0“]2s 
laminate.  In  the  case  of  the  laminate,  this  trend  was 

opposite  to  the  ones  proposed  by  the  three  two-parameter  models. 

Notched  Strength  Predictions 

The  two  two-parameter  models  of  Nulsmer  and  Whitney  utilize  the 
exact  anisotropic  infinite  plate  solution  for  centrally  located  cracks 
and  holes.  As  a  result,  comparisons  of  these  models  with  the  inherent 
flaw  ittdel  and  BIE  method  can  only  be  made  for  center  notched  (CN) 
laminates.  Comparisons  between  experimental  and  predicted  strengths  of 
the  CN  specimens  for  different  aspect  ratios  and  laminate  orientations 
are  as  shown  in  Figures  B-12  to  B-15  for  our  T300/934  graphite/epoxy 
materials.  Equations  {B-4).  (B-6)  and  (B-7)  of  the  two-parameter 
nodels  were  used  to  obtain  the  analytical  curves  representative  of 
those  iodels.  The  characteristic  lengths  used  to  generate  these 
curves  were  assumed  to  be  constant  and  specific  values  as  suggested  by 
the  oHginators  of  the  models  were  used  as  was  discussed  earlier. 

As  may  be  noted  in  Figures  B-12  to  B-15,  considerable  differences  were 
found  between  predictions  based  upon  the  three  two-parameter  models 
and  the  experimental  results.  While  some  of  the  variations  may  have 
been  due  to  experimental  error,  it  is  felt  that  the  fundamental  reason 
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Fiqure  B-13.  Comparison  Between  Experimental  and  Analytical  Critical 
Remote  Stresses  for  CN  [0°/90  Laminates. 
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for  the  difference  beween  the  various  methods  of  approach  were 
likely  due  to  the  assumed  constant  characteristic  lengths.  Obviously, 
tetter  comparison  between  analysis  and  experiment  would  have  been  ob- 
Ulned  If  a  -best  fit"  characteristic  length  had  been  used.  Such  a 
procedure  is  nothing  more  than  "curve  fitting"  and  does  not  lead  to 
a  better  understanding  of  the  fracture  behavior  of  coctiosites. 

The  results  labelled  Cruse  in  Figures  B-12  to  B-15  were  obtained 
by  assuming  that  the  fracture  toughness  calculated  using  the  aspect 
ratio.  2a/w  •  0.25.  represented  the  true  value  for  the  laminate. 

This  value  was  chosen  because  in  all  cases  the  fracture  toughness 
approached  asymptotic  values  beyond  2a/w  •  0.25.  I.  these  figures, 
the  best  Ciwarison  between  analysis  and  experiment  was  seen  using  the 

one-parameter  BIE  model  of  Snyder  and  Cruse. 

The  three  two-parameter  fracture  models  can  be  used  for  circular 
holes  as  well  as  for  crack-like  flaws.  Thus,  using  the  sane  constant 
characteristic  lengths  as  discussed  above  for  the  AS-3501  graphite/ 
epoxy  specimens  conUining  holes,  comparisons  between  the  predicted 
fracture  strengths  of  the  three  two-parameter  models  and  experiments 
are  shown  in  Figures  B-16  and  B-17  for  a  variety  of  laminate  orienta¬ 
tions.  Generally,  the  inherent  flaw  predictions  correlated  better 
with  the  experimental  fracture  strengths  than  the  two  approaches  of 
Mhitney  and  Nuismer.  Again,  a  better  correlation  between  analysis  and 
eiiperiment  could  have  been  obtained  by  adjusting  the  characteristic 
dimension  which  again  would  have  represented  a  "curve  fitting"  process 
This  was  not  done  for  the  same  reasons  discussed  earlier. 


CRITICAL  GROSS  STRESS,  o.  (ksl) 


A  [A5°]q5 
□  t90“]85 

WHITNEY,  ET  AL. 
-  POINT  STRESS 

-  AVERAGE  STRESS 

-  WADDOUPS,  ET  AL 


10^  § 


Figure  B-16.  Comparison  of  Experimental  and  Theoretical  Fracture 
figure  Strengths  for  AS-3501  Graphite/Epoxy  Laminates 

Containing  Centrally  Located  Holes.  (Results 
corrected  to  infinite  width  using  isotropic  stress 
concentration  factor.) 
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Figure  B-17. 
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Finite  Width  Correction  Factor 

When  correlating  experimental  results  uttl.  analyses  based  on  a 
body  «1th  infinite  dimensions  or  yice-.ersa.  finite  width  correction 
factors  (FHCF)  must  be  employed  if  the  correlations  are  to  be  meaning¬ 
ful.  To  date,  most  of  the  FUCF's  for  isotropic  materials  are 
MaKrically  obtained.  In  the  case  of  composite  materials,  all  the 
FWCfs  are  obuined  numerically.  Thus,  a  comparison  of  the  isotropic 
niCF  and  the  anisotropic  FHCF  (obtained  from  the  BIE  method)  for  the 
nouhed  specimens  used  in  this  investigation  would  be  appropriate. 

Since  the  BIE  method  is  a  numerical  method  and  like  other 
„«erical  methods,  its  convergence  should  be  determined.  The  procedure 
adopted  here  was  to  detemine  the  segment  spacings  first,  followed  by 
prescribing  the  boundary  conditions.  The  results  obuined  were 
a.,«red  to  the  results  reported  1n  reference  [BO]  for  isotropic 
maurlals.  For  the  three  types  of  implanted  notches  used  here,  the 
difference  between  the  BIE  results  and  the  ones  in  reference  [BO]  was 
less  than  0.3t.  Upon  completion  of  this  preliminary  work,  the 
conversion  to  anisotropic  cases  was  just  a  matter  of  defining  the 
laBinate  orientation  and  principal  material  properties. 

The  variations  between  the  Isotropic  FHCF  (equation  (B-ll))  and 
the  anisotropic  values  for  the  CN  specimens  are  as  shown  in  Figure 
B-18.  Except  for  the  [tdS'],^  laminate,  the  difference  between  the 
isotropic  and  anisotropic  values  are  less  than  5*  for  all  aspect  ratios. 
In  the  case  of  the  [i45»]4s  laminate,  the  most  anisotropic  laminate 
investigated  here,  the  maximum  difference  between  the  Isotropic  and 
anisotropic  value  was  approximately  10*  when  2a/w  =  0.45. 
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Comparisons  between  the  isotropic  FWCF  [80]  and  the  anisotropic 
FWCF  of  the  DEN  and  SEN  laminates  are  as  shown  in  Figures  B-19  and 
8-20.  respectively.  Variations  in  the  FWCF's  were  noted,  especially 
in  the  DEN  cases.  For  these  two  notches,  the  largest  variations 
between  the  isotropic  and  anisotropic  FWCF's  were  seen  for  the 
laminate.  Quantitatively,  the  maximum  deviation  was  about  20%  at  the 

largest  aspect  ratio. 

Thus,  depending  on  the  anisotropy,  aspect  ratio  and  type  of  im¬ 
planted  notch,  the  isotropic  FWCF  may  or  may  not  vary  substantially 
from  its  anisotropic  counterpart.  However,  it  should  be  noted  that 
the  isotropic  FWCF's,  for  the  three  types  of  implanted  notches  in¬ 
vestigated  here,  were  much  easier  and  convenient  to  compute  than  the 

anisotropic  FWCF's. 

Conclusion  and  Future  Considerations 
From  the  results  and  discussion  here  on  the  fracture  behavior  of 
graphite  epoxy  laminates,  the  following  observations  can  be  made: 

.  The  COO  measurements  of  the  CN  specimen  indicated  that 

the  compliance  of  square  slots  was  always  higher  than  narrow 
slits,  for  all  laminates  and  aspect  ratios.  Consequently, 
the  remote  fracture  stresses  of  square  slots  were  the  higher 

of  the  two  cases. 

.  [±45*]4s  laminate  containing  circular  holes  had  a  slightly 

lower  remote  fracture  stress  than  those  containing  narrow 
slits  even  though  the  compliance  of  the  holes  was  higher. 
This  conclusion  was  drawn  from  observations  of  the 


F(a/vi) 


photoelastic  coatings  which  revealed  that  the  intense 
regions  were  larger  in  the  case  cf  the  holes.  Thus,  holes 
are  more  critical  than  slits  for  this  laminate  geometry. 

The  effects  of  anisotropy  on  the  stress  field  in  the  vicinity 
of  the  notch  tips  of  the  CN  specimens  were  clearly  illustrated 
by  the  birefringent  coatings.  By  varying  the  anisotropy  of 
the  material  from  quasi -isotropic  to  highly  anisotmoic. 
the  pattern  of  the  fringes  was  varied.  For  laminates, 

the  effects  of  interlamina  constraints  were  detected. 

Stable  crack  growth  was  observed  in  most  of  the  in-situ 
microscopic  observations  and  was  generally  associated  with 
matrix-dominated  situations.  Microcracks  were  always  present 
In  the  notched  and  unnotched  free  surfaces  of  the  failed 
general  laminates* 

For  CN  laminates,  when  non-self-similar  crack  propagation 
occurred,  characteristic  lengths  decreased  as  the  aspect  ratio 
Increased.  When  self-similar  crack  propagation  occurred,  the 
characteristic  lengths  approached  an  asymptotic  value  as  the 
aspect  ratio  increased,  which  was  also  seen  in  the  results 
of  the  work  of  the  originators  of  the  two-paraneter  models. 
Using  the  suggested  characteristic  dimension  of  0.04  inch 
(1.016  m)  for  the  point  stress  and  inherent  flaw  models  and 
0.15  inch  (3.810  mm)  for  the  average  stress  model,  a 
comparison  of  the  predicted  fracture  strength  with  the  experi¬ 
mental  results  of  specimens  containing  CN  and  holes  was  made. 
For  CN  specimens  the  best  correlation  between  analysis  and 
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experiment  was  obtained  using  the  BIE  method.  The  inherent 
Haw  model  generally  correlated  belter  with  experimental 
results  than  the  other  models  used  for  circular  hole 
geometries.  No  "curve-fitting"  was  attempted  with  the  two- 

parameter  models. 

.  The  finite  width  correction  factors  used  to  correlate  the 
experimental  results  with  analysis  varied  according  to  the 
anisotropy  of  the  laminate.  The  maximum  deviation  between 
the  isotropic  and  anisotropic  correction  factors  (obtained 
from  the  BIE  method)  was  approximately  lOX  for  the 
CK  laminate  at  2a/w  »  0.45  and  approximately  20?  for  the 
£45*]^^  SEN  and  DEN  cases  when  a/w  *  0.7  and  2a/w  ■  0.7, 
respectively. 

In  conclusion,  the  BIE  method  is  an  excellent  analytical  tool  for 
investigating  the  fracture  behavior  of  composite  materials.  The  two- 
parameter  Bodels  and  isotropic  finite  correction  factors  should  only 
be  used  as  a  rule  of  thumb  type  judgment  or  design  situation  and 

prinrily  only  for  quasi -isotropic  laminates. 

Tentatively,  future  investigations  on  the  time-dependent  fracture 
behavior  of  notched  composites  should  utilize  the  rationale  of  the  BIE 
Kthod.  This  is  primarily  because  it  is  more  convenient  to  use  and 
has  a  better  accuracy  than  other  numerical  methods  such  as  finite 
elements  and  finite  difference.  The  results  and  discussion  presented 
here  should  be  sufficient  justification  for  the  preliminary  time- 
dependent  analysis  of  notched  composites  using  the  BIE  method. 
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THE  TIME-TEMPERATURE  BEHAVIOR  OF  GRAPHITE/ EPOXY  LAMINATES 

by 

Yew  Thye  Yeow 

(ABSTRACT) 

An  in-depth  study  on  the  time- tempera tore  behavior  of  epoxy 
based  continuous  and  elastic  fiber  reinforced  composite  materials  and 
the  feasibility  of  using  the  time-temperature  analogy  as  an  accelerated 
characterization  method  to  predict  the  long-term  behavior  are  presented. 
This  is  a  two-pronged  investigation  as  the  material  investigated 
(graphite/epoxy)  is  essentially  quasi-elastic  at  room  temperature  and 
viscoelastic  at  elevated  temperatures.  Correlations  between  analysis 
and  experiments  are  presented  whenever  possible. 

At  elevated  temperatures,  master  curves  of  matrix  dominated  uni¬ 
directional  laminates  are  obtained  using  the  time-temperature  super¬ 
position  principle.  Using  the  principal  properties  (E,^,  £22(0. 

G^^Ct)  and  v^2^»  "^^ter  curves  of  other  off-axis  laminates  reduced  at 
any  arbitrary  ambient  temperature  are  predicted.  The  results  obtained 
from  the  short-tem  (16-minute)  tests,  medium  tenn  (25-hour)  tests  and 
predictions  are  shown  to  correlate  reasonably  well.  In  addition,  the 
delayed  failure  prediction  and  experimental  results  are  shown  to 
correlate  reasonably  well.  Both  the  experimental  and  predicted  delayed 
failure  results  indicate  that  the  creep  strength  master  curves  are 
dependent  on  the  in-plane  stress  states  and  that  failure  modes  vary 


fmn  ductile  to  ouusi-brittle  depeodin,  on  the  duretion  of  the  tests. 

At  room  temperature,  a  sumnary  of  the  unnotched  and  notthed 
behavior  is  presented.  In  the  unnotched  case,  total  stress-strain 
responses  of  general  syimnetric  laminates  are  predicted  «ith  a  non- 
linear  analysis  and  compared  with  experimental  results.  The  analysis 
ts  used  to  validate  intralamina  shear  test  methods.  In  the  notched 
case,  test  results  tor  specimens  containing  variations  in  notch 
^tries  and  anisotropy  are  given.  These  results  are  correlated  to 
three  two-parameter  and  one  one-parameter  analytical  fracture  models. 


